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ABSTRACT 

Context. Young loose nearby associations are unique samples of close (<150pc), young («5-100Myr) pre-main sequence (PMS) 
stars. A significant number of members of these associations have been identified in the SACY collaboration. We can use the proxim¬ 
ity and youth of these members to investigate key ingredients in star formation processes, such as multiplicity. 

Aims. With the final goal to better understand multiplicity properties at different evolutionary stages of PMS stars, we present the 
statistics of identified multiple systems from 113 confirmed SACY members. We derive multiplicity frequencies, mass-ratio, and 
physical separation distributions in a consistent parameter space, and compare our results to other PMS populations and the field. 
Methods. We have obtained adaptive-optics assisted near-infrared observations with NACO (ESO/VLT) and IRCAL (Lick 
Observatory) for at least one epoch of all 113 SACY members. We have identified multiple systems using co-moving proper-motion 
analysis for targets with multi-epoch data, and using contamination estimates in terms of mass-ratio and physical separation for targets 
with single-epoch data. We have explored ranges in projected separation and mass-ratio of a [3-1000 au], and q [0.1-1], respectively. 
Results. We have identified 31 multiple systems (28 binaries and 3 triples). We derive a multiplicity frequency (ME) of 
MF3_ioooau=28.4]^3g% and a triple frequency (TF) of TF3_ioooau=2.8]^gg% in the separation range of 3-1000au. We do not find any 
evidence for an increase in the ME with primary mass. The estimated mass-ratio of our statistical sample (with power-law index 
7=-0.04 ± 0.14) is consistent with a flat distribution (y = 0). 

Conclusions. Analysis from previous work using tight binaries indicated that the underlying multiple system distribution of the SACY 
dataset and the young star-forming region (SFR) Taurus are statistically similar, supporting the idea that these two populations formed 
in a similar way. In this work, we show further similarities (but also hints of discrepancies) between the two populations: flat mass- 
ratio distributions and statistically similar ME and TF values. We also compared the SACY sample to the field (in the separation range 
of 19-100 au), finding that the two distributions are indistinguishable, suggesting a similar formation mechanism. 
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1. Introduction 

The multiplicity of a stellar population holds a huge amount 
of information about the birth, evolution and eventual fate of 
the population itself and its individual members. Thus, compre¬ 
hensively characterising these properties has implications across 
many fields such as star formation, N-body dynamics, and su¬ 
pernovae rates. Multiplicity surveys of solar-type stars (e.g. 
Duquennoy & Mayor 1991; Raghavan et al. 2010; Tokovinin 
2014) have made use of a combination of techniques (from small 

* Based on observations obtained using the instruments NACO at 
the VET (077.C-0483, 081.C-0825, 088.C-0506, 089.C-0207) and the 
adaptive optics facilities at the Lick observatory (UCO 3m) 

** Appendix A is available in electronic form at http: //www. aanda. 
org. Table B. 1 is only available at the CDS via anonymous ftp to http: 
//cdsarc.u-strasbg.fr(ftp://123.45.678.9) or via http:// 
cdsarc.u-strasbg.fr/viz-bin/qcat?3/A+A/XXX/XXX 


physical separations to large: spectroscopy, interferometry in¬ 
cluding sparse aperture masking, - SAM -, speckle imaging, 
adaptive optics, - AO -, assisted imaging, and classical direct 
imaging) to probe and characterise multiple systems with peri¬ 
ods between ~1-10* ** ° days (with a peak at :«10^ day). As a re¬ 
sult, they show that almost one out of every two solar-type star 
(0.9-1.5 Mq) has one or more companion. Considering the high- 
mass end of the initial mass function (IMF), the O-type stars, the 
percentage increases up to ^100% (Sana et al. 2014). In other 
words, these massive stars form almost exclusively in multiple 
systems. 

Such surveys are volume-limited, making use of the high 
number of main sequence (MS) field stars that are nearby to pro¬ 
duce robust statistics (e.g., in the case of Tokovinin 2014, 90% 
completeness, considering their spatial density, for d < 70 pc). 
However, the history of this population is not well constrained 
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as it comprises stars bom at different times and in a range of 
environments. As a result, they must have undergone quite dif¬ 
ferent dynamical processing(Parker et al. 2009). Therefore, the 
information regarding the multiplicity properties as a direct out¬ 
put from star formation has been smeared out during its some¬ 
what utrreproducible dynamical history. To study a more pris¬ 
tine population one has to consider the observations of younger 
stars (<100 Myr) in their nurseries, that usually provide lower 
number statistics. Such studies are usually performed in nearby 
star-forming regions (SFRs) like e.g. Taurus (Leinert et al. 1993; 
Kraus et al. 2011; Daemgen et al. 2015), Ophiuchus/L1688 
(Ratzka et al. 2005), Chamaeleon I (Lafreniere et al. 2008), 
IC 348; (Duchene 1999), or the ONC (Reipurth et al. 2007). 
However, as King et al. (2012b) demonstrated, the amount of 
overlapping cluster-to-cluster parameter space limits the num¬ 
ber of systems that one is able to statistically compare (typically 
< 100 ). 

Young associations have also been studied due to their prox¬ 
imity (<150 pc) and youth (<100Myr). There are a handful of 
studies focused on the stellar multiplicity in these populations 
(see e.g. Brandeker et al. 2003, 2006). The results of these 
works lead to varying conclusions about the nature of the stars 
in such associations. As an example, Brandeker et al. (2003) 
found a very high frequency of stellar multiple systems in TW- 
Hydrae, and an absence of systems (with separations >20 au) 
in rj Chamaeleontis, with an extremely small likelihood (< 10 "^) 
that the two populations share the same parental multiplicity dis¬ 
tribution. It is still unclear whether this is due to dependencies 
on purely initial density (Moraux et al. 2007) or additional envi¬ 
ronmental effects (Becker et al. 2013). 

There have also been some studies focused on the substel- 
lar multiplicity in these associations (Masciadri et al. 2005; 
Lafreniere et al. 2007; Nielsen & Close 2010; Evans et al. 2012; 
Chauvin et al. 2010; Biller et al. 2013). In addition to such sur¬ 
veys, many individual objects have been imaged with the aim to 
disentangle the planetary / brown dwarf nature of their compan¬ 
ions and further understand their circumstellar discs (Chauvin 
et al. 2004; Huelamo et al. 201 1). 

However, a detailed statistical analysis of the abundance and 
characteristics of multiple systems in young associations com¬ 
pared to those identified in SFRs is still missing. Such an anal¬ 
ysis was very complex until recently, posing a compromise be¬ 
tween the publicly available data and the ability to probe the 
same parameter space. Motivated by this, in 2006 we started an 
homogeneous observational campaign to derive the multiplicity 
of young loose associations. 

We aim to characterise the multiplicity properties of these 
associations across a large parameter space {a ~ 0.1-10"^ au) by 
combining a range of observational techniques. We have already 
derived the spectroscopic multiple system fraction (the tight¬ 
est multiple systems, semi-major axis (a) <1 au) of our sample 
(Elliott et al. 2014). We found that the derived multiplicity frac¬ 
tion is compatible with that of the field and nearby SFRs, as one 
would expect due to the lack of dynamical processing for sys¬ 
tems with such high binding energy, that have statistically sim¬ 
ilar primordial abundances. We are now extending our analysis 
to wider binaries, using AO imaging. 

AO imaging allows us to identify and characterise multiple 
systems across a huge parameter space for targets that belong 
to nearby associations. As an example, for an average member 
of our sample (c/=70pc, age=30Myr, Mi-1 Mq), an instrument 
such as NAOS-CONICA at the VET (point spread function, - 
PSF-, full-width half maximum, -FWHM-, of 0.07 " at 2.2 fim, 
with a strehl ratio of 1), is sensitive to companions at separa¬ 


tions [a~3-1000au] in the [M 2 ~ 0.01-1 Mq] mass range. As 
mentioned previously, the peak of field stars period distribution 
is at 10^ day which, considering a secondary of 0.1 Mq, corre¬ 
sponds to a semi-major axis of ~40 au assuming a circular orbit. 
Therefore, our observing campaign should yield a high number 
of multiple systems when a large enough sample is observed. 
In addition, the observations must have sufficient temporal cov¬ 
erage to discard that a visual pair is just a chance alignment 
through a proper motion analysis (as orbital motion is negligi¬ 
ble in the majority of cases). 

The manuscript is organised as follows; Section 2.1 de¬ 
scribes the target selection. Section 3 details the observations, 
data reduction, and the sensitivity limits, as well as the literature 
analysis. Section 4 outlines the methods used to identify bound 
multiple systems and the uncertainty estimations. Section 5 de¬ 
scribes how we derived the multiplicity properties of our sample. 
Section 6 presents the statistical analysis of the built distributions 
and comparisons to other populations. Section 7 states the main 
conclusions of this work and prospects. 


2. Sample definition 

The search for associations containing young stars (SACY) 
dataset contains 2000 stars located all over the sky. The stars 
are defined as members of distinct moving groups by the conver¬ 
gence method (see census of Torres et al. 2006). The kinematic 
properties of the stars are derived from the parallax, the proper 
motion (PM) and the radial velocity (RV) values. 

The starting point for our sample selection was based on the 
census by Torres et al. (2006). For the southern targets two ad¬ 
ditional criteria were followed; (i) we only consider stars with 
spectral types later than G5, in order to increase the contrast be¬ 
tween the primary and the possible companion candidate, and 
(ii) the right ascension (RA) should be obtainable from Paranal 
observatory when our service observations were conducted (pe¬ 
riods 77, 81, 88, and 89). The northern targets also followed cri¬ 
terion (i) and were also observed according to their visibility 
during the two runs at the Lick observatory in 2008 and 2009. 

Later on, the membership of several targets was revised 
based on new data (e.g. new RV values derived from high- 
resolution spectra, Elliott et al. (2014). This revision allowed 
us to eliminate contamination within the sample from interlop¬ 
ers. From the 201 SACY targets observed originally with AO, 
113 are now confirmed SACY members, while 88 have been re¬ 
jected. 

The bulk of this work is focused on the 113 confirmed SACY 
members, and the specific details of the sample are provided be¬ 
low (summarized in Table 1). However, for completeness with 
our original observations, we also provide details on the 88 non- 
SACY targets, including their observational parameters and a 
summary of the identified multiple systems, in Tables B.l and 
C.l, respectively. 

2.1. Filtered coherent sample for this work 

Table 1 provides the census of the 113 confirmed SACY mem¬ 
bers (once the additional data from Elliott et al. 2014 have been 
included to identify interlopers) studied in this work. We have 
included two K magnitude values; the one from the 2MASS sur¬ 
vey, K (2MASS), and our estimated K (Mi). From our obser¬ 
vations (with a narrow field of view that does not allow for a 
robust flux calibration), we can only obtain relative photome¬ 
try and therefore we use 2MASS magnitudes for calibration. We 
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denoted our estimated magnitude, K (Mi), and it only differs 
from the original 2MASS value in the case of a detection of a 
new companion (previously unresolved) within the 2MASS PSF. 
The typical PSF width for 2MASS is ai2"and, in some cases, 
multiple components can be found within this angular distance. 
To estimate the photometry of the primary and components we 
weighted the 2MASS value according to the observed flux ratio. 
In the next sections, the primary masses have been calculated 
using this value, and the age value included in Table 2 with the 
evolutionary tracks of Baraffe et al. (1998). 

Finally, we note that all the targets in this sample have 
PM values from the UCAC4 catalogue (Zacharias et al. 2012). 
Precise PMs will allow us to determine whether identified com¬ 
panions are co-moving with the primary source. 
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Table 1. All SACY targets observed and analysis in this work. The distance values are kinematic, derived from the convergence method (Torres 
et al. 2006). 


ID 

RA 

hh;mm:ss.s 

DEC 

dd:mm:ss 

K (2MASS) 
(mag.) 

K(Mi) 

(mag.) 

Ass. 

Distance 

(pc) 

Ml 

(Mo) 

CD-45 14955B 

23:11:53.6 

-45:08:00 

8.85 

8.85 

ABD 

49.3 

0.56 

GJ 9809 

23:06:04.8 

-h63:55:34 

6.98 

6.98 

ABD 

24.5 

0.61 

HD 1405 

00:18:20.9 

-h30:57:22 

6.51 

6.51 

ABD 

27.3 

0.78 

HD 160934 

17:38:39.6 

-h61:14:16 

6.81 

6.81 

ABD 

33.1 

0.8 

HD 201919 

21:13:05.2 

-17:29:12 

7.58 

7.58 

ABD 

39.5 

0.72 

HD 222575 

23:41:54.2 

-35:58:39 

7.62 

7.62 

ABD 

63.6 

0.96 

HD 24681 

03:55:20.4 

-01:43:45 

7.25 

7.25 

ABD 

53.3 

0.96 

LP 745-70 

16:33:41.6 

-09:33:11 

7.55 

7.55 

ABD 

30.3 

0.59 

TYC 486-4943-1 

19:33:03.7 

-h03:45:39 

8.66 

8.66 

ABD 

72.0 

0.77 

V* PX Vir 

13:03:49.6 

-05:09:42 

5.51 

5.51 

ABD 

21.7 

0.9 

HD 152555 

16:54:08.1 

-04:20:24 

6.36 

6.36 

ABD 

46.8 

1.13 

GJ 885 A 

23:00:27.9 

-26:18:42 

6.27 

6.75 

ABD 

32.7 

0.81 

Wolf 1225 

22:23:29.1 

-h32:27:34 

6.05 

6.77 

ABD 

15.1 

0.48 

BD-03 4778 

20:04:49.3 

-02:39:20 

7.92 

8.06 

ABD 

71.5 

0.91 

GJ4231 

21:52:10.4 

-h05:37:35 

7.38 

8.2 

ABD 

30.6 

0.5 

TYC 91-82-1 

04:37:51.4 

-h05:03:08 

8.65 

8.86 

ABD 

87.6 

0.82 

CD-52 9381 

20:07:23.7 

-51:47:27 

7.39 

7.39 

ARG 

29.4 

0.59 

2MASS J05200029H-0613036 

05:20:00.2 

-h06: 13:03 

8.57 

8.57 

BPC 

67.8 

0.67 

BD-13 6424 

23:32:30.8 

-12:15:51 

6.57 

6.57 

BPC 

27.9 

0.7 

CD-31 16041 

18:50:44.4 

-31:47:47 

7.46 

7.46 

BPC 

53.2 

0.87 

CD-54 7336 

17:29:55.0 

-54:15:48 

7.36 

7.36 

BPC 

68.4 

1.06 

CPD-72 2713 

22:42:48.9 

-71:42:21 

6.89 

6.89 

BPC 

37.3 

0.81 

GSC 07396-00759 

18:14:22.0 

-32:46:10 

8.54 

8.54 

BPC 

95.2 

0.93 

HD 161460 

17:48:33.7 

-53:06:43 

6.78 

6.78 

BPC 

71.0 

1.2 

Smethells 20 

18:46:52.5 

-62:10:36 

7.85 

7.85 

BPC 

52.4 

0.73 

TYC 6872-1011-1 

18:58:04.1 

-29:53:04 

8.02 

8.02 

BPC 

82.6 

0.99 

V* TX PsA 

22:45:00.0 

-33:15:25 

7.79 

7.79 

BPC 

20.1 

0.19 

V* V4046 Sgr 

18:14:10.4 

-32:47:34 

7.25 

7.25 

BPC 

76.9 

1.15 

V* WW PsA 

22:44:57.9 

-33:15:01 

6.93 

6.93 

BPC 

20.1 

0.35 

GJ 3322 

05:01:58.7 

-h09:58:59 

6.37 

6.8 

BPC 

37.8 

0.85 

CD-26 13904 

19:11:44.6 

-26:04:08 

7.37 

7.77 

BPC 

78.9 

1.03 

CD-27 11535 

17:15:03.6 

-27:49:39 

7.38 

8.06 

BPC 

87.3 

1.01 

GSC 08350-01924 

17:29:20.6 

-50:14:53 

7.99 

8.65 

BPC 

66.3 

0.63 

2MASS J05203182H-0616115 

05:20:31.8 

-h06:16:11 

8.57 

8.7 

BPC 

71.0 

0.66 

2MASS J08110934-5555563 

08:11:09.3 

-55:55:56 

9.4 

9.4 

CAR 

129.0 

0.84 

2MASS J08521921-6004443 

08:52:19.2 

-60:04:44 

9.37 

9.37 

CAR 

163.9 

1.04 

2MASS J08563149-5700406("> 

08:56:31.4 

-57:00:40 

9.76 

9.76 

CAR 

191.6 

0.94 

BD-03 5579 

23:09:37.1 

-02:25:55 

7.82 

7.82 

CAR 

63.0 

0.85 

CD-44 1533 

04:22:45.6 

-44:32:51 

8.58 

8.58 

CAR 

103.7 

0.91 

CD-49 4008 

08:57:52.1 

-49:41:50 

8.64 

8.64 

CAR 

102.9 

0.89 

CD-53 2515 

08:51:56.4 

-53:55:56 

8.75 

8.75 

CAR 

137.1 

1.13 

CD-54 2499 

08:59:28.7 

-54:46:49 

8.4 

8.4 

CAR 

109.5 

1.08 

CD-55 2543 

09:09:29.3 

-55:38:27 

8.4 

8.4 

CAR 

115.2 

1.12 

CD-61 2010 

08:42:00.4 

-62:18:26 

8.83 

8.83 

CAR 

121.4 

0.93 

CD-69 783 

10:41:23.0 

-69:40:43 

8.38 

8.38 

CAR 

86.5 

0.87 

HD 107722 

12:23:29.0 

-77:40:51 

7.14 

7.14 

CAR 

59.0 

1.05 

HD 309751 

09:31:44.7 

-65:14:52 

8.36 

8.36 

CAR 

141.6 

1.2 

HD 8813 

01:23:25.8 

-76:36:42 

6.79 

6.79 

CAR 

46.5 

0.92 

TYC 8174-1586-1 

09:11:15.8 

-50:14:14 

9.5 

9.5 

CAR 

117.2 

0.79 

TYC 8557-1251-1 

07:55:31.6 

-54:36:50 

9.19 

9.19 

CAR 

123.0 

0.86 

TYC 8582-3040-1^“^ 

08:57:45.6 

-54:08:36 

9.35 

9.35 

CAR 

143.2 

0.89 

TYC 8962-1747-1 

11:08:07.9 

-63:41:47 

8.29 

8.29 

CAR 

88.2 

0.9 

TYC 9217-417-1 

09:59:57.6 

-72:21:47 

8.69 

8.69 

CAR 

83.6 

0.8 

TYC 9486-927-1 

21:25:27.4 

-81:38:27 

7.34 

7.34 

CAR 

36.3 

0.71 

HD 22213 

03:34:16.3 

-12:04:07 

6.79 

6.97 

CAR 

52.4 

0.94 

BD-07 2388 

08:13:50.9 

-07:38:24 

6.92 

7.4 

CAR 

42.6 

0.77 

2MASS J09131689-5529032 

09:13:16.8 

-55:29:03 

8.36 

8.61 

CAR 

131.7 

1.15 

CD-57 1709 

07:21:23.7 

-57:20:37 

8.7 

8.7 

CAR 

100.5 

0.87 
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ID 

RA 

hh mm ss 

DEC 
dd mm ss 

K (2MASS) 
(mag.) 

K(Mi) 

(mag.) 

Ass. 

Distance 

(pc) 

Ml 

(Mo) 

CPD-62 1293 

09:43:08.8 

-63:13:04 

8.6 

8.82 

CAR 

69.4 

0.7 

2MASS J08371096-5518105 

08:37:10.9 

-55:18:10 

9.38 

9.38 

CAR 

180.8 

1.11 

2MASS J03573723-0416159 

03:57:37.2 

-04:16:15 

8.75 

8.75 

COL 

106.5 

0.89 

2MASS J09322609-5237396 

09:32:26.0 

-52:37:39 

8.84 

8.84 

COL 

101.9 

0.85 

BD-H08 742 

04:42:32.0 

-h09:06:00 

9.12 

9.12 

COL 

105.6 

0.81 

BD-11 648^'=) 

03:21:49.6 

-10:52:17 

9.26 

9.26 

COL 

128.1 

0.86 

CD-36 1785 

04:34:50.7 

-35:47:21 

8.59 

8.59 

COL 

80.1 

0.8 

HD 26980 

04:14:22.5 

-38:19:01 

7.62 

7.62 

COL 

80.7 

1.12 

HD 27679 

04:21:10.3 

-24:32:20 

7.81 

7.81 

COL 

78.0 

0.94 

HD 31242 

04:51:53.5 

-46:47:13 

8.16 

8.16 

COL 

68.1 

0.82 

TYC 9178-284-1 

06:55:25.1 

-68:06:21 

8.94 

8.94 

COL 

108.3 

0.86 

V* V479 Car 

09:23:34.9 

-61:11:35 

7.96 

7.96 

COL 

87.7 

1.07 

2MASS J02303239-4342232 

02:30:32.4 

-43:42:23 

7.23 

7.57 

COL 

51.3 

0.81 

V* V1221 Tau 

03:28:14.9 

+04:09:47 

7.44 

8.01 

COL 

81.5 

0.92 

2MASS J08240598-6334024 

08:24:05.9 

-63:34:02 

8.13 

8.18 

COL 

118.6 

1.2 

HD 272836 

04:53:05.1 

-48:44:38 

8.24 

8.34 

COL 

78.6 

0.84 

BD-16 351 

02:01:35.6 

-16:10:00 

7.96 

8.43 

COL 

80.9 

0.83 

2MASS J04272050-4420393 

04:27:20.4 

-44:20:39 

8.56 

8.7 

COL 

86.2 

0.81 

CD-43 1395 

04:21:48.6 

-43:17:32 

8.42 

9.16 

COL 

141.0 

0.93 

2MASS J11594226-7601260 

11:59:42.2 

-76:01:26 

8.3 

8.3 

ECH 

101.0 

0.98 

2MASS J12194369-7403572 

12:19:43.6 

-74:03:57 

8.86 

8.86 

ECH 

108.9 

0.78 

2MASS J12392124-7502391 

12:39:21.2 

-75:02:39 

7.78 

7.78 

ECH 

108.6 

1.2 

CD-69 1055 

12:58:25.5 

-70:28:49 

7.55 

7.55 

ECH 

106.7 

1.2 

HD 104467 

12:01:39.1 

-78:59:16 

6.85 

6.85 

ECH 

106.5 

1.2 

V* DZ Cha 

11:49:31.8 

-78:51:01 

8.49 

8.49 

ECH 

102.7 

0.91 

V* MP Mus 

13:22:07.5 

-69:38:12 

7.29 

7.29 

ECH 

106.1 

1.2 

HD 105923 

12:11:38.1 

-71:10:36 

7.18 

7.26 

ECH 

117.1 

1.2 

TYC 9245-535-1 

12:56:08.3 

-69:26:53 

7.99 

8.34 

ECH 

119.3 

1.14 

2MASS J12202177-7407393 

12:20:21.7 

-74:07:39 

8.37 

8.4 

ECH 

114.8 

1.07 

2MASS J05581182-3500496 

05:58:11.8 

-35:00:49 

9.38 

9.38 

OCT 

105.8 

0.77 

2MASS J06400573-3033089 

06:40:05.7 

-30:33:08 

8.7 

8.7 

OCT 

162.4 

1.2 

CD-58 860 

04:11:55.6 

-58:01:47 

8.36 

8.36 

OCT 

87.1 

0.88 

CD-66 395 

06:25:12.3 

-66:29:10 

9.0 

9.0 

OCT 

126.1 

0.91 

HD 271037^“) 

05:06:50.5 

-72:21:11 

8.67 

8.67 

OCT 

145.8 

1.19 

BD-18 4452A 

17:13:11.6 

-18:34:25 

6.48 

6.48 

OCT 

16.7 

0.51 

CD-30 3394A 

06:40:04.9 

-30:33:03 

8.59 

8.59 

OCT 

162.4 

1.2 

CD-47 1999 

05:43:32.1 

-47:41:10 

8.64 

8.68 

OCT 

166.7 

1.2 

BD-20 nil 

05:32:29.3 

-20:43:33 

8.7 

8.7 

OCT 

129.7 

1.1 

HD 274576 

05:28:51.3 

-46:28:19 

8.81 

8.92 

OCT 

117.1 

0.89 

2MASS J04302731-4248466 

04:30:27.3 

-42:48:46 

8.73 

9.42 

OCT 

122.8 

0.82 

2MASS J06033540-4911256 

06:03:35.4 

-49:11:25 

9.08 

9.81 

OCT 

173.8 

0.89 

2MASS J05490656-2733556 

05:49:06.5 

-27:33:55 

8.25 

8.25 

THA 

74.3 

0.86 

ASAS J051536-0930.8 

05:15:36.4 

-09:30:51 

8.08 

8.08 

THA 

77.5 

0.91 

BD-12 943 

04:36:47.1 

-12:09:20 

7.76 

7.76 

THA 

62.3 

0.88 

BD-19 1062 

04:59:32.0 

-19:17:41 

8.07 

8.07 

THA 

64.1 

0.83 

CD-35 1167 

03:19:08.6 

-35:07:00 

7.72 

7.72 

THA 

44.6 

0.72 

CD-46 1064 

03:30:49.0 

-45:55:57 

7.1 

7.1 

THA 

42.2 

0.84 

CD-53 544 

02:41:46.8 

-52:59:52 

6.76 

6.76 

THA 

41.4 

0.9 

CD-78 24 

00:42:20.3 

-77:47:39 

7.53 

7.53 

THA 

50.1 

0.83 

TYC 8098-414-1 

05:33:25.5 

-51:17:13 

8.16 

8.16 

THA 

51.6 

0.68 

TYC 9344-293-1 

23:26:10.7 

-73:23:49 

7.94 

7.94 

THA 

44.2 

0.65 

HD 22705 

03:36:53.4 

-49:57:28 

6.14 

6.18 

THA 

43.3 

1.04 

CD-44 1173 

03:31:55.6 

-43:59:13 

7.47 

7.58 

THA 

44.0 

0.75 

TYC 8083-455-1 

04:48:00.6 

-50:41:25 

7.92 

8.25 

THA 

54.8 

0.7 

2MASS J05182904-3001321 

05:18:29.0 

-30:01:32 

8.3 

8.41 

THA 

66.2 

0.76 


Notes. Targets excluded from statistical analysis due to the quality of their sensitivity curves, see Figure 5. 
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3. Observations and data reduction 

3.1. Southern Targets: NACO data 

We had four sets of observations in service mode dedicated 
to collecting data on SACY sources with NAOS-CONICA, an 
adaptive optics system and near-lR camera at the very large tele¬ 
scope (VLT). The first epoch of data was taken between March- 
August 2006 (077.C-0483), and 76 targets were observed (see¬ 
ing ssO.9")- The 40 objects (out of the 76) identified to have 
potential companions were re-observed (seeing ssO.7") between 
March-June 2008 (081.C-0825). We had two additional pro¬ 
grams to observe other members of the SACY associations in 
2011 and 2012 (088.C-0506, and 089.C-0207), and we observed 
60 and 25 targets (seeing =sl.0"), respectively (see Table A.l for 
details on individual sources). 

All the targets are bright enough (V < 16 mag or K < 13 mag) 
to be used as AO reference stars. For the majority of the stars, we 
used the visible wave-front sensor (VIS WFS). For the reddest 
objects we used the IR WFS with the N90C10 dichroic, which 
sends 90% of the light to NAOS and 10% to the infra-red detec¬ 
tor. Since all the sources are very bright in the near-IR, the ob¬ 
servations with the VIS WFS were performed with intermediate- 
and narrow-band filters to avoid saturation of the primary star. 
The first selection of filters was IB 2.21 (do, A/l=2.21, 0.06//m) 
and NB2.17 (do, Ad=2.17, 0.023/rm). During the observations, 
the IB 2.21 filter showed strong optical ghosts that prevented us 
from detecting companions in some particular regions of the de¬ 
tector. Therefore, we replaced it with the IB 2.27 (do, Ad=2.27, 
0.06 jum) filter. The observations with the IR WFS were per¬ 
formed with the broad band Kj filter (do, Ad=2.18, 0.35 ;um). 
All the observations from 2006 and 2008 used the S27 objec¬ 
tive, which provides a field of view (FOV) of 27.6" x 27.6". The 
data from 2011 and 2012 used the S13 objective, with a FOV of 
14"X14". 

We observed each object using a jitter sequence to remove 
the sky contribution and to correct for bad pixels. The sky emis¬ 
sion was estimated from the median image of all offsets posi¬ 
tions. The average total exposure time for each source was =sl2 
minutes. We reduced the data using the recommended eclipse 
ESO software (Devillard 1997) and the esorex driven interface 
gasgano. Dark subtraction, bad pixel correction, and flat-held 
correction were applied to each frame prior to aligning and 
stacking. 

To derive accurate astrometry of the potential companions, 
we derived the plate scale and orientation on the infrared detec¬ 
tor, using archival observations of the astrometric field 6^ Ori 
C (McCaughrean & Stauffer 1994), and the astrometric binaries 
IDS 2150 and IDS 22141 (van Dessel & Sinachopoulos 1993); 
see Chauvin et al. (2012) for details. The values measured in dif¬ 
ferent epochs are displayed in Table 3. For each target, we have 
adopted the plate scale and orientation closest to the date of ob¬ 
servation. 

A summary of the NACO observations is provided in 
Table 4. 

3.2. Northern Targets: Lick data 

To observe the 40 northern targets in our sample, we used the 
adaptive optics facility at the Lick Observatory Shane 3 m tele¬ 
scope (Gavel et al. 2002) between 14-18 July 2008 and 8-9 July 
2009. 

In 2008 the conditions were relatively good, although not 
photometric, with an average seeing between 0.8" and 1.1" in 


Table 3. Detector plate scales (i + CTj) and true north orientations 

(t ± o-t) 


Obs. date 

Plate scale 

True north" 

Target 


(mas / pix) 

(deg) 


NACO 

24/04/2006 

27.01±0.05 

-0.08±0.14 

e' Ori C 

26/08/2006 

27.02±0.05 

-0.12±0.20 

e' Ori C 

06/04/2008 

27.01±0.05 

-0.36±0.14 

e' Ori C 

10/06/2008 

27.01±0.05 

-0.39±0.15 

IDS 2150, 




0' Ori C 

15/11/2011 

13.27±0.05 

-0.47±0.14 

0* Ori C 

04/07/2012 

13.27±0.05 

-0.78±0.20 

IDS 22141 

Lick 

15/07/2008 

75.7±0.1 

-tl.80±0.30 

HIP 66225 

09/07/2009 

75.7±0.2 

-0.30±0.20 

HIP 66225 


Notes. This value is subtracted from the measured position angle 
(PA). 


the visible, at the zenith. The targets were all bright enough to 
be used as a reference for the wave-front sensor. All of the tar¬ 
gets were observed using the Kj filter, which provides the best 
compromise between the AO performances and the sensitivity to 
faint and red companions. 

When a clear companion was detected in the images, an ad¬ 
ditional H band (and sometimes J band) image was acquired to 
provide colour information. All observations were performed us¬ 
ing a dither pattern to efficiently remove cosmic rays, bad pixels 
and accurately compute the sky. The data were processed using 
Eclipse, including dark subtraction, flatfield correction, sky sub¬ 
traction, frame registration, and stacking. A Hipparcos binary 
(HIP 66225, Gontcharov 2006) was observed every night to be 
used as an astrometric calibrator and derive an accurate plate 
scale and plate position angle on the sky (see Table 3). 

The 2009 observations aimed at obtaining second epoch im¬ 
ages of the 21 targets with companion candidates detected in 
the first epoch images. Conditions were poorer, with a seeing 
of 1.5" and some cirrus. As we were most interested in de¬ 
riving accurate astrometry, the observations were performed in 
a single band (using the K, filter). The same observing strategy 
(five point dither pattern) and processing steps were applied. The 
same astrometric reference was used to calibrate the plate scale 
and position angle of the camera (see Table 3). 

A description of the Lick data used in this work is included 
in Table 4. 

3.3. Archival data and literature search 

We queried the ESO archive for related NACO data for objects 
that were observed only once in our service campaigns. This 
way, we could utilise a second epoch to determine if the de¬ 
tected companion candidates are co-moving. In total 64 objects 
had multi-epoch observations. In addition to this, we queried the 
CDS for object-type determinations and multiplicity flags for all 
201 targets studied in this work, (see Table A.l). Eor objects pre¬ 
viously identified as multiple systems, we checked whether the 
companions’ properties were within the parameter space of our 
observations, which in turn, tests the completeness of our work. 
As a result, we did not find any previously identified systems in 
this parameter space which we failed to detect. 
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Table 2. Summary of the properties and frequency of visual multiple systems identified for the SACY associations studied in this 
work. The minimum and maximum ages represent the dispersion found in the literature. 


Name 

ID 

Dist. 

(pc) 

Age*"' 

(Myr) 

Max. / Min. 
(Myr) 

N. of observed**' 
objects 

N. ofVBs*"' 
(this work) 

MFs-ioooau*'^* 

(%) 

AB Doradus 

ABD 

34±26 

100 

100 / 50 

16(16) 

5 

31/^^ 

Argus 

ARG 

I06±51 

40 

50/20 

1 

0 

- 

/3-Pic MG 

BPC 

31±2I 

22 

30/ 10 

17 (17) 

5 

29/» 

Carina 

CAR 

85±35 

35 

35/10 

26 (24) 

5 

2Lf 

Columba 

COL 

82±30 

35 

40/ 15 

17(16) 

6 

38/» 

e-Cha 

ECH 

I08±9 

9 

9/6 

10(10) 

2 

20/^’ 

Octans 

OCT 

I41±34 

38 

38/10 

11(10) 

5 

50/ij 

Tuc-Hor 

THA 

48±7 

33 

33/10 

14(14) 

3 

21/^'' 


Notes. Ages are derived from the convergence method described in Torres et al. (2006) using new radial velocity values, results to appear in 
Torres et al. in prep.**' Effective number used to calculate ME values shown in parentheses.*'^' Systems either confirmed from PM analysis or high 
probability (>95% detections) based on single-epoch data. *“*' In the separation range 3-1000 au considering observed targets with sensitivities 
>95% of our sample. 


Table 4. A summary of all AO-imaging data used in this work. 
The number of observations refers to the number of reduced im¬ 
ages as opposed to separate objects, i.e., an object can be counted 
twice if it was re-observed. 


Instrument 

FOV 

(arcsec) 

Filter 

^0 

(//m) 

AT 

(pm) 

No. of 
Obs. 

NACO 

14 X 14 

Ks 

2.18 

0.35 

91 


28 X 28 

Ks 

2.18 

0.35 

53 


28 X 28 

IB 2.27 

2.27 

0.06 

82 


28 X 28 

IB 2.21 

2.21 

0.06 

11 


28 X 28 

NB 2.17 

2.17 

0.02 

12 

Lick 

20x20 

Ks 

2.15 

0.32 

70 


20x20 

H 

1.66 

0.30 

14 


20x20 

J 

1.24 

0.27 

2 


3.4. Observational detection limits 

Figure 1 shows the average contrast in the selected biters as a 
function of the angular separation from the primary star for both 
the NACO (using the S13 and S27 objectives) and Lick data. On 
average, we reach AKj = 7 mag and AKs = 4 mag at an angular 
separation of 0.5", for the NACO and Lick data, respectively. For 
an age of 30Myr (the average age of the SACY associations) 
and a primary mass of LOMq, these differences in magnitude 
correspond to masses of ~15Mjup and ~0.10Mq (dashed lines 
in Figure 1), using the evolutionary tracks from Baraffe et al. 
(1998). 

4. Identifying multiple systems 

In order to identify multiple systems in our sample, we have 
followed the next procedure; brst, we have identibed all the 
point sources within the FOV of each target. Then, we have de¬ 
termined the astrometry and photometry of these sources, and 
used these parameters to determine whether the identibed point 
sources are bound companions or background objects. We pro¬ 
vide all the details in the next subsections. 



Fig. 1. The average sensitivity limit in all biters for both NACO 
and Lick data. The limits were computed from the 3cr root-mean 
squared (RMS) of the PSF radial proble for each target, and then 
averaged within each biter. The Ks biter for the NACO data is 
displayed for the S13 and S27 objective. The dashed lines rep¬ 
resent our mass limits at (X'S, according to Baraffe et al. (1998) 
evolutionary tracks, assuming an average age of 30Myr and a 
primary mass of 1 Mq. 


4.1. Initial source detection 

Combining the NACO data and Lick data, we have imaged a 
total of 201 young targets. In order to identify companions in 
our images, we brst need to identify all possible point sources 
in the FOV. We used the source extraction software SExtractor 
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(Bertin & Amouts 1996) to perform the source identification. 
Only detections of sources with 5 or more pixels with values 
> 3cr were considered (where cr is the standard deviation of the 
background signal). 

We visually checked each image to (i) conhrm whether the 
detection was from a true point source as opposed to, for exam¬ 
ple, an optical ghost and (ii) search for potential severely blended 
binaries (with component separations ~4 pixel) that had been de¬ 
tected as one source. With the catalogues of point sources for 
each observation we then modelled the PSFs using the software 
PSFEx in conjunction with SExtractor. Eor the majority of ob¬ 
servations, the primary targets were used to model the PSE. In the 
case that the primary target was resolved into a close binary, ei¬ 
ther further components in the EOV (as close to the centre of the 
image as possible) were used, or another bright target observed 
during the same night with the same instrumental setup and as 
close in time as possible. Eigure 2 shows a typical example of a 
blended object originally detected as a single source. The cen¬ 
ter and the photometric aperture as derived by SExtractor are 
represented by the black circle and black ellipse, respectively. 
Through the modelling of the PSE we were able to extract the 
astrometry and photometry for the two components, whose cen¬ 
tres and EWHMs are shown by the red crosses and dashed cir¬ 
cles, respectively. 

101 -,- 



Ql -.- 

0 5 10 

Pixel 

Fig. 2. A section of the NACO image for target CD-30 3394B, 
showing two partially resolved components. The black solid- 
lined ellipse is the aperature applied by SExtractor and the 
black filled circle the centre of the ellipse. The red crosses and 
black dashed circles indicate the centre and EWHM of the model 
PSE extracted from the primary CD-30 3394A, using PSEEx for 
each source. 


4.2. Relative photometry 

Our observations did not include photometric calibrator helds 
and, therefore, we can only perform relative photometry in our 
images. However, the vast majority of our targets were imaged 
in the Ks band, so we can make use of the available 2MASS 
Ks values (Skrutskie et al. 2006). In the cases where we discov¬ 


ered multiple components within the original PSE of 2MASS, 
we weighted the magnitudes according to the flux ratios of the 
components. These values are provided in the hfth column of 
Table 1, as K(Mi), and have been used to compute the mass ratio 
of the multiple systems using evolutionary tracks (Baraffe et al. 
1998), and the age estimation from Table 2. 

4.3. Relative astrometry 

For every point source detected in each observation, we calcu¬ 
lated the angular separation from the primary in detector coordi¬ 
nates and then applied the appropriate pixel scale and true north 
correction value (Table 3). In the case of non-blended sources, 
we have used the parameters from the aperture extraction, while 
for severely blended sources we have used the ones coming from 
PSF modelling. We have performed both types of extractions 
(aperture and PSF) on a small number of targets and confirmed 
that, for non-blended sources, the resultant photometry and as¬ 
trometry is consistent within the uncertainties. We then have 
employed a different method to determine whether the detected 
point sources are bound companions or not, depending on the 
availability of multiple epochs. 

4.3.1. Uncertainties 

The main source of uncertainty comes from the temporal and 
spatial variability of the atmospheric conditions. The uncertain¬ 
ties for the angular separation, p(x,y) + (Tp, and position angle, 
6(x, y) ± cTg, were derived using equations 1 and 2, using the ap¬ 
propriate plate-scale and true north correction values. Table B.l 
shows the parameters and their uncertainties of the multiple sys¬ 
tems identified for each observation. 



where ctx = xcts, o-y - yo-g and ctj is the uncertainty in the 
plate scale. 



where cr, is the uncertainty in the true north correction. 
4.3.2. Multi-epoch data 

To assess whether a pair of point sources are bound, and not a 
projection effect, we computed their relative motion and com¬ 
pared it to the motion one would expect for a background ob¬ 
ject. The PM values were taken from UCAC4 (Zacharias et al. 
2012). We then combined this with the parallactic motion, and 
produced the total relative motion one would expect to see if the 
point source were a background, stationary object (see Chauvin 
et al. 2005). An example is shown in Figure 3. We would expect 
a bound companion to show negligible movement with respect 
to the primary between our multi-epoch observations, and this 
is our criterion for binary identification. In those cases in which 
the PM of the source did not produce significant motion, given 
the time difference between observations, we used the criterion 
described in Section 4.3.3. 


8 










P. Elliott et aL: Search for associations containing young stars (SACY) 


140 

120 


100 


80 

60 
to 
<1 

40 

20 

0 

-20 

-40 -30 -20 -10 0 10 20 30 40 

A a (mas) 

Fig. 3. The relative movement of four point sources, in RA and 
DEC, for two observations of GSC 08350-01924. The black dot¬ 
ted line and black markers represent the motion one would ex¬ 
pect from a background (BKG) source with respect to the pri¬ 
mary, and the shaded area represents the error space. The dark 
blue marker represents a source showing no significant relative 
motion, consistent with a bound (BND) companion at 0.3". The 
other coloured markers represent three background sources de¬ 
tected at separations of 10.5, 13.2 and 14.2" from the primary. 



4.3.3. Single-epoch data 

For those targets with only one epoch of data (NACO data using 
S13 camera) we could use their properties, to provide statistical 
constraints on how likely is the source to be bound. One way 
to do that is to estimate the potential contamination from back¬ 
ground sources given the target’s galactic co-ordinates, and the 
limiting magnitude of the observations. However, galactic mod¬ 
els such as Robin et al. (2003) have a limited resolution and, 
given our small FOV (13"x 13"), this method only provides a 
very rough estimate of the level of contamination. To avoid this 
limitation and provide more realistic estimations, we decided to 
use our available multi-epoch data to assist in the classification 
of our single-epoch data. 

We classified our own multi-epoch detections as either 
bound companions or background sources from PM analysis, 
see Section 4.3.2). We then used two additional sources of data. 
Firstly, Wahhaj et al. (2013), a high-contrast direct imaging sur¬ 
vey for giant planets. We converted the H-band photometry into 
mass-ratios using the evolutionary tracks of Baraffe et al. (2003) 
and the ages quoted in the work. Secondly, we used the data from 
Daemgen et al. (2015), a multiplicity survey in Taurus conducted 
in Ks-band. From these two works, we only used sources classi¬ 
fied by PM analysis (in total, 246 classified sources). 

We then pseudo-randomly split this classified data: 2/3 as 
a training set, 1/3 a test set. Using support vector machines 
(SVM)', we determined the optimum soft-boundary in the phys¬ 
ical separation-mass-ratio space that splits bound and back- 


* see Chapter 9 of Ivezic et al. (2014) for many applications of this 
technique in classification of astronomical objects such as differentiat¬ 
ing variable and non-variable MS-stars using photometric colours. 


ground sources. On average, we reached a completeness and 
contamination rate of >80% and <10% respectively, when clas¬ 
sifying our test set. With this classifier we then calculated the 
probability that our single-epoch sources are bound (Pbnd) or 
background (Pbkg)- In this instance, there are only two labels: 
either bound or background and as a result, Pbnd = 1 - Pbkg- We 
repeated this procedure 1000 times and took the median result to 
account for the variations induced by different training and test 
datasets. The individual probabilities (Pbnd) of these companion 
candidates with their respective properties are shown in Table 5. 

To test whether this technique produces realistic probabili¬ 
ties, we calculated the bound probabilities for the 35 companions 
already confirmed by PM analysis. The median bound probabil¬ 
ity was ss99%. 

Figure 4 shows the graphical result of this procedure. 
There is a clear trend towards high mass-ratio, closely-separated 
sources being bound (as one would expect), and only very few 
bound sources lie in the parameter space mainly populated by 
background sources. 



0.2 0.4 0.6 0.8 1 

Mass-ratio 


Fig. 4. The probability of identifying a bound system from 100- 
0% in intervals of 10% from blue to red in physical separation- 
mass-ratio parameter space. The circles and crosses represent the 
training and test set, respectively. Blue markers represent bound 
companions, and red markers unbound companions. 


5. Results 

We have identified 31 potential multiple systems (28 binaries 
and 3 triples) from 113 confirmed SACY targets. Of these 31 
systems, 7 have been confirmed by PM analysis and 24 have a 
bound-probability >95%^ (see Section 4.3.3). The details of the 
observational properties from individual observations are dis¬ 
played in Table B.l. A summary of the physical properties of 
the multiple systems is included in Table 5. 

We have followed the multiple component designations de¬ 
scribed in Tokovinin (2005), where systems are described by 
individual components and super components. For example, a 
triple system composed of a binary system orbiting the primary 
(brightest) star would have the following designations: ”A, B, 
and ”Ba, Bb, B” following increasing hierarchy. In such cases, 
there are two entries in Table 5, one to describe ”A, Ba” and the 
other ”Ba, Bb”. 

The analysis presented in the following subsections consid¬ 
ers the 31 identified multiple systems from 8 SACY associations 

^ Based on this statistical analysis there is a 57% probability that all 
24 systems are bound from the product sum of their respective proba¬ 
bilities. 
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Fig. 5. Top panel: The physical projected separations versus 
mass-ratio for identified multiple systems. The filled red mark¬ 
ers represent binary companions, the other coloured markers are 
components of triple systems with matching colours for the sec¬ 
ondary and tertiary. The size of the marker is proportional the 
companion’s mass. The grey lines represent the sensitivity limits 
for the individual targets. The red and black lines represent the 
90% and 95% detections limits for the whole sample, respec¬ 
tively. Bottom Panel: The 95% detection limits for four studies, 
red: this work, blue: Daemgen et al. (2015), green: Kraus et al. 
(2011) and brown: Raghavan et al. (2010), used for analysis in 
Section 6 . 


(outlined in Table 2). We note that have also identified 28 mul¬ 
tiple systems in the sample of non-SACY members. A summary 
of these systems’ properties can be found in Appendix C. 

We have defined two ranges of physical separations in our 
analysis, derived from the sample sensitivities displayed in 
Figure 5. The first, 3-1000 au, is the most extensive, consider¬ 
ing 109 from the 113 observed targets, and including all the 34 
companion detections (28 binaries and 3 triples). We define an¬ 
other as 10-1000 au, that will allow us to compare our results 
with other studies. For this parameter space, we again consider 
109 observed targets but only 27 companions detections (25 bi¬ 
naries and 1 triple). 

5 .1. Multiplicity frequencies 

The frequency of multiple systems in a population can be com¬ 
puted in a number of ways, see Reipurth & Zinnecker (1993). 
In this paper we consider the following quantities: multiplicity 
frequency (MF), companion star frequency (CSF) and triple fre¬ 
quency (TF). 

We have identified 31 multiple systems from observations of 
113 confirmed SACY targets. Therefore, we estimate a raw mul¬ 
tiplicity frequency (MFja„) of 27.4^^^%, and a raw companion 
star frequency (CSFraw) of 30. 1^3 ^%. However, as mentioned 
previously, these frequencies include companions that would not 
be detected around all of the targets due to different contrasts 
and distances (see Figure 5). When we remove these targets, we 
derive a multiplicity frequency of MF3_ioooau=28.4;|;3 g% in the 
separation range 3-1000 au. All of our companion detections are 
above 95% sensitivity curve and therefore we do not disregard 
any of them in our calculation. 


5.1.1. Higher-order multiple systems 

Considering the higher-order systems separately (3 triple sys¬ 
tems) we derive a triple frequency of TF 3 _ioooau= 2 . 8 ;^Qg% in 
the separation range of 3-1000 au. It is interesting to note that 
2 of the 3 triple systems identified in this work have secondary 
components at very small separations (4.8 and 5.2 au), that are 
only detectable due to the small distance to the systems (31 and 
17 pc). 

5.1.2. Multiplicity frequency and primary mass 

We analysed the effect of the primary mass on the multiplic¬ 
ity frequency within our sample for the 109 systems in the sep¬ 
aration range 3-1000 au. We binned our data according to the 
Freedman-Diaconis rule, which is suited for non-Gaussian dis¬ 
tributions (bin width^O. 1 Mq). The results are shown in Figure 6 . 
The large uncertainties for masses of 0.25, 0.35 and 0.45 Mq are 
a result of the low-number fractions, 0 / 1 , 0/1 and 1 / 1 , respec¬ 
tively. 



Fig. 6. Primary mass versus multiplicity frequency for systems 
in the separation range 3-1000 au. The data is binned according 
to the Freedman-Diaconis rule. The fractions are shown by each 
mass bin. 


5.2. Mass-ratios 

We calculated the mass-ratios (q) of our multiple systems in two 
ways, to test the effect of using primary mass values derived 
with different methods. We note that, in the case of triple 
systems composed of a primary component A and an outer 
orbiting binary (Ba, Bb), the mass-ratio of the component Ba is 
derived as be Mea/MA: 

(i) We derived primary mass values (Mi) from absolute 
magnitudes using kinematic distances (D) and ages with the 
evolutionary tracks of Baraffe et al. (1998). We then calculated 
the companion masses from AK values and computed the 
mass-ratios (the normalised cumulative distribution function 
(CDF) is shown by the red line in Figure 7). 

(ii) We used the AK values of our components for a range 
of Ki values, which is equivalent to a range of model-based 
Mi values for the given age of the system (0.8-1.2Mo, accord¬ 
ing to the evolutionary tracks of Baraffe et al. 1998). To ob¬ 
tain the uncertainty, we calculated the standard deviation of the 
mass-ratio as a function of the primary mass for a given sys¬ 
tem. We then computed 100 mass-ratio distributions composed 
of pseudo-random realisations of each mass-ratio value and its 
uncertainty. We binned the data in 0.01 steps and took the av- 
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erage and standard deviation producing a smoothed, normalised 
CDF, shown by the black dots and shaded areas (dark and light; 
Icr and 3cr) of Figure 7, respectively. 

The values and uncertainties used in further analysis (shown 
in Table 5) are obtained using method (ii) in order provide real¬ 
istic uncertainties given that, at this time, the masses are model- 
based. 


Fig. 7. Left panel. Normalised frequency of mass-ratios. The 
dashed line represents a power-law function with index y - 0. 
Right panel: The cumulative frequency distribution of mass- 
ratios using the two methods outlined in Section 5.2. Method 
(i) The red line represent the derived mass ratios. Method (ii) 
The black dots, dark and light blue-coloured area represent the 
average, 1 cr and 3 cr variation in frequency from 100 realisations 
of the mass ratio, respectively. 


The mass distribution is usually described using a power law 
index, y. We derive a value of y = -0.04 + 0.14 using linear 
regression of 1000 pseudo-random realisations of the mass-ratio 
distribution considering their respective uncertainties. 

5.3. Physical separations 

We calculated the physical separations of the multiple systems 
using the kinematic distances shown in Table 1, and the av¬ 
erage angular separation from our observations. The top panel 
of Figure 5 shows the physical separation versus mass-ratio for 
the 31 SACY multiple systems, considering both multiple-epoch 
and single-epoch observations with probabilities >95%). The 
size of the markers represents the companion’s mass. In addi¬ 
tion to this, the detection limits are shown for each multiple sys¬ 
tem (grey lines) and 90% and 95% detection limits for the entire 
observed sample in red and black, respectively. This calculation 
demonstrates the powerful nature of the SACY dataset when ob¬ 
served with AO-imaging; we can probe very small physical pro¬ 
jected separations across our sample down to low mass-ratios. 

In addition to this, we generated separation distributions 
in two ranges. Firstly, 10-1000 au, for SACY, Taurus (Kraus 
et al. 2011; Daemgen et al. 2015, hereafter Kll and D15) and 
the field (Raghavan et al. 2010, hereafter RIO), see Figure 8. 
Secondly, 19-100au, for SACY, Taurus (only Kll) and the field 
(see Figure 9). The first separation range was chosen to max¬ 
imise the number of comparable systems due to different obser¬ 
vational techniques. The second, to make the same comparison 
that was made in King et al. (2012a). 

First of all, we created 90% sensitivity curves for each study 
in terms of mass-ratio and physical separation (shown by the 
coloured lines in the bottom panel of Figure 5), using the avail¬ 
able published data. In the case of RIO, we used the sensitivity 
curve from Figure 11 of that paper. Then, we defined the mas¬ 


ter sensitivity curve as the upper boundary of all the individual 
curves (shown by the black solid line in the bottom panel of 
Figure 5). We then only consider systems with parameters above 
this curve. 

To perform this analysis we used kernel density estimations 
(KDEs), see Scott (2009) for details. This technique minimises 
potential problems with the bin size and the phase. The results 
for the two separation ranges are shown in the left panels of 
Figures 8 and 9. The plots shows the KDE from SACY in red, 
Taurus in blue and green, and the field in gold with 95% confi¬ 
dence intervals (CIs) shown by the respective-coloured shaded 
areas. The CIs were calculated by considering the distribution 
of densities for each physical separation from 10^ bootstrap it¬ 
erations of the KDE. To see the differences between the distri¬ 
butions more clearly we have plotted them in the form of cu¬ 
mulative frequency distributions (CDEs) in the right panels of 
Eigures 8 and 9. 

The low-density nature of Taurus makes it very interesting to 
compare to our associations, as they most likely formed in a low- 
density environment too. We use both the work of K11 and D15 
since both samples are complementary. The statistical sample 
defined in D15 consists of 10 members from the extended Taurus 
region (~20 Myr) and 5 from the young Taurus region (~2 Myr), 
totalling 15 multiple systems. That of Kll only considers the 
young region and in our defined parameter space, this produces 
a sample of 50 multiple systems for comparison. 

The field represents a far more processed population of stars 
from a range of initial environments and, therefore, provides a 
mixed sample of evolved multiple systems. The number of com¬ 
parable systems in the field within our defined parameter range 
is 123. 



Fig. 8. Left: Red, green, blue and gold shaded areas repre¬ 
sent the 68% confidence intervals for the kernel density esti¬ 
mation (KDE) for SACY, Taurus (Kll), Taurus (D15) and the 
field (RIO) respectively. The multiple systems considered are all 
above the master sensitivity curve defined in the bottom panel 
of Eigure 5 and in the separation range 10-1000 au. Right: The 
cumulative distribution frequencies for the same systems. 


Eirstly, we tested the effect of primary mass on the separa¬ 
tion distribution as the studies we are considering probe slightly 
different primary mass ranges. 

We performed a set of simulations using 10,000 synthetic bi¬ 
naries with primary masses between 0.2-1.4Mo. We performed 
two sets of simulations which differed by their mass-ratio distri¬ 
bution. In the first case, we considered the mass-ratio distribution 
as uniform across all primary masses. In the second simulation 
we adopt the information from Eigure 5 of K11 as priors. In the 
range Mp 0.7-1.4Mo the mass-ratio distribution as uniform. Eor 
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the range Mi: 0.2-0.7Mo there is a preference for more equal 
mass-ratio systems for lower primary masses. We created a prob¬ 
ability function for the secondary masses, and pseudo-randomly 
chose mass-ratio values from this distribution. 

The period values were drawn from the log-normal distri¬ 
bution dehned in Raghavan et al. (2010), ju=5.03 and cr-2.28 
log (day). We then htted the separation distributions with log¬ 
normal distribution and evaluated their parameters. When com¬ 
paring the two sets of simulations, we see that the peaks of the 
separation distributions are shifted to smaller separations for the 
lightest primary masses: mean values of 30 au in the hrst case, 
and 60 au in the second case. However, this difference is not 
significant given the separation range we are considering ( 10 - 
1000 au) There is also no signihcant difference resulting from 
the mass-ratio distributions. 

Based on this result, we have compared separations distribu¬ 
tions of held stars in the mass range of 0.7-1.4Mo, and SACY 
objects in the mass range of 0.2-1.2 Mq. 



Fig. 9. Left: Red, green and gold shaded areas represent the 68 % 
conhdence intervals for the kernel density estimation (KDE) of 
multiple systems for SACY, Taurus (Kll) and the held (RIO), 
respectively. The multiple systems considered are all above the 
master sensitivity curve dehned in the bottom panel of Figure 5 
and in the separation range 19-100 au. Right: The cumulative 
distribution frequencies for the same systems. 


Our comparison between the three populations is shown in 
Figure 8 . The shape of the SACY distribution differs from those 
of the other populations at a 1 cr level ( 68 % Cls are represented 
by the coloured hlled areas in left panel of Figure 8 ). The distri¬ 
bution departs in form from the other populations at aiSO-lOO au. 
Potentially, this departure for wider systems is the result of a bias 
in the way we classify our objects as bound, favouring close-in 
systems, using the method described in Section 4.3.3. However, 
to check this, we included companions regardless of their proba¬ 
bilities and the resultant distribution still differed from the other 
populations at a 1 cr level for the widest systems. 

6. Discussion 

The SACY sample offers a unique opportunity to study multi¬ 
plicity among young stars down to very small physical separa¬ 
tions. Currently, it is not known whether all these associations 
shared a common origin, due to limitations on the accuracy of 
their galactic motion^. Whether these associations share a com- 

^ The GAIA mission (Lindegren et al. 2008) will provide extremely 
precise astrometry for members of such associations, with expected ac¬ 
curacies of 7-25 pas. 


mon origin or not, it is likely that these groups resulted from 
sparse star formation whereby fewer than 100 stars were formed 
in each group. Therefore by studying these targets we are obtain¬ 
ing crucial information regarding this mode of star formation. In 
this analysis we present the results assuming a common forma¬ 
tion mechanism, i.e. grouping the data together in some cases 
from all associations and looking for relationships in physical 
parameters (i.e. primary mass, physical separation). 

6.1. Frequency of multiple systems 

The triple frequency is TF 3 _ioooau= 2 . 8 j^Qg% in the separation 
range 3-1000 au. If we consider the range 10-1000 au, we derive 
the value TFio-ioooau= 0 . 9 );Q 3 %. The value from Daemgen et al. 
(2015) in Taurus, TFio_ioooau=l-81j'5%, is comparable within 
the uncertainties. 

The multiplicity frequency for SACY is 
MF3_ioooau=28.4;|;3g% in the separation range 3-1000au 
across the primary mass range 0.2-1.2Mo. To compare to the 
work of Daemgen et al. (2015) in the Taurus SFR, we consider 
the separation range 10-1000 au and derive a multiplicity fre¬ 
quency of MFio-ioooau=25.7j;3y%. The multiplicity frequency 
for Taurus in that parameter space is MFio-ioooau=26.3);4g%, 
so the two values agree very well within the uncertainties. 
However, if we limit the primary mass of our targets the result 
is slightly different, as discussed below. 

As shown by King et al. (2012b,a), collating studies from 
Taurus, p-Ophiuchus, Chamaeleon I, IC 348 and the ONC, and 
additionally the work of Kraus et al. (2011); Daemgen et al. 
(2015) studying Taurus, the multiplicity frequency of Taurus 
is marginally higher (^1 cr) than that of other SFRs. In order 
to make a comparison between Taurus (Daemgen et al. 2015), 
SACY and the held (Raghavan et al. 2010) we have selected tar¬ 
gets in the separation range 10-1000 au and the primary mass 
range >0.7 Mq. In the case of SACY the primary mass range is 
0.7-1.2 Mq for Taurus and the held it is 0.7-1.4 Mq. The resul¬ 
tant multiplicity frequencies (MFio_ioooau,Mi>o. 7 ) are shown in 
the bottom panel of Figure 10. The result of SACY (23.2^3 3 %) 
agrees with the held (28.9)^2 9 *^) and marginally disagrees with 
Taurus (37.9!®-^%) at =sl cr level. Based on these results, within 
this range, there is marginal evidence that the Taurus region has 
a higher frequency of multiple systems then SACY and the held. 
This result further supports the idea that Taurus has a higher 
MF value than other SFRs and associations. Kroupa & Bouvier 
(2003) showed how environments such as Taurus are insuffi¬ 
cient at binary disruption and, therefore, their high MF values 
are due to their binary populations being more pristine than that 
of heavily processed populations. We would therefore expect, as¬ 
suming similar primordial multiple system distributions, to see 
the same effect in the SACY associations. Currently, due to low 
numbers, we cannot say conclusively whether this property is 
shared or not. There is no disagreement between SACY and 
the held, which, assuming common primordial multiple distri¬ 
butions, suggests that few multiple systems have been destroyed 
by dynamical processing in this range. 

In contrast to the general trend of high multiplicity in low- 
density SFRs and associations, Brandeker et al. (2006) found 
an absence of multiple systems with separations >20 au in the 
young ( 6-8 Myr, Jilinski et al. 2005) cluster 77 Chamaeleontis. 
This survey of 17 targets used NACO observations similar to 
those presented here, probing similar targets and mass-ratio sys¬ 
tems. They found a probability <10 that 77 Chamaeleontis and 
TW-Hydrae share the same parent multiple system distribution. 
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Fig. 10. Age versus multiplicity frequency. Top panel: In the 
separation and primary mass range 3-1000 au and 0.2-1.2Mo, 
respectively. Bottom panel: In the separation and primary mass 
range range 10-1000 au and 0.7-1.4Mo, respectively. Data from 
(Daemgen et al. 2015, Taurus) and (Raghavan et al. 2010, the 
field) are shown by the star and diamond marker, respectively. 
Argus (ARG) has been omitted as only 1 object was observed. 


Although 77 Chamaeleontis is denser than the associations pre¬ 
sented in this work, the authors conducted analysis to disregard 
dynamical processing as the cause of such a result, finding a 
collision time-scale 3000 times that of the cluster age. The fre¬ 
quency of spectroscopic systems (SBs) in tj Chamaeleontis is 
statistically similar to that of the SACY associations (Riviere 
et al. submitted), =sl0%. The physical separations of these SBs 
are so tight that they should be free of dynamical processing 
and therefore similarities between regions support the idea of 
statistically similar binary populations. However, the discs in rj 
Chamaeleontis are long-lived (Sicilia-Aguilar et al. 2009), con¬ 
sidering its age, with respect to other young regions. Therefore, 
at this time, it is still not clear whether the differences in mul¬ 
tiplicity frequency between this cluster and other low-density 
SFRs and associations, are the result of dynamical processing 
or stem from environment-dependent star formation. 

However, since there is a general trend of high multiplic¬ 
ity frequency with low-density environments, it is likely that, at 
least in some part, it stems from the lack of dynamical process¬ 
ing compared to stars in populations such as the field (Lada & 
Lada 2003). 

6.1.1. What is the effect of the primary mass on the 
frequency? 

It is common in many populations to observe an increase in 
multiplicity frequency with primary mass (see Figure 12 of 
Raghavan et al. 2010). We also noted the same relationship 
among the tightest multiple systems, i.e. spectroscopic bina¬ 
ries (Elliott et al. 2014). However, this is considering a very 
large mass range, from spectral types O-M. We have investigated 
whether this relationship is also observable within our sample of 
stars (primary mass range 0.2-L2Mo.) 


Figure 6 shows the result of binning our data in terms of 
primary mass for the 109 targets with >95% sensitivity from 
Figure 5). There is no observed trend of increasing multiplicity 
frequency with primary mass. There is a slight disagreement be¬ 
tween some mass bins (^1 cr) but only in the case of very low 
numbers (< 3). The most reliable mass bins (with more than 
10 objects) show no tendency of higher frequency with primary 
mass. This is in agreement with the results of Kraus et al. (201 1) 
and Daemgen et al. (2015) within the same primary mass range. 

6.2. The impact of a flat mass-ratio distribution 

Duchene & Kraus (2013) compiled observations and discussed 
how random pairing from the IMF (Chabrier 2003) does not 
agree with the derived observational mass-ratio distributions for 
primary masses of ~l-10Mo. The mass-ratio power-law index 
is a much steeper function of the mass for lower-mass primaries 
(< 0.5 Mo). 

The mass-ratio distribution derived for the SACY sample 
(power-law index 'y=-0.04 + 0.14), agrees with the compilation 
of observational results in Duchene & Kraus (2013), in the pri¬ 
mary mass range of =s0.2-L2Mo (see left panel of Figure 7). 
The derived y values in this range are very close to 0, i.e. flat 
distributions. With our current observations we are not sensitive 
to more extreme mass-ratio systems (smaller than 0 . 2 ) at close- 
separations. To investigate whether the power-law index contin¬ 
ues to be flat in this primary mass range, we need to observe our 
sample with more sensitive techniques, such as extreme high- 
contrast imaging or sparse-aperture masking (Evans et al. 2012). 
We aim to perform this analysis in future work. 

6.3. Is there dynamical evolution within the SACY sample? 

Another great advantage of the SACY sample is the range of 
PMS ages that we can probe (~10-100Myr). This allows us to 
search for any potential dynamical evolution within the sample 
on these time scales. 

Figure 10 shows the multiplicity frequency for systems in 
the separation range 3-1000 au. If we consider the youngest and 
oldest associations, e-Cha (ECH) and AB Dor (ABD), we do not 
see significant differences in the multiplicity frequency. Given 
the number of observations we have at this time, we are limited 
by low-number statistics. To investigate this further, we need to 
observe more systems in each association to build up larger sam¬ 
ples. Given that there is only one association that does not agree 
with the other seven, we assume a statistically similar primordial 
multiple system distribution for all associations. 

Although the parameters of such systems may evolve on a 
time-scale of 10-100 Myr (see below) it is very unlikely that the 
systems will be entirely destroyed. Most of the destruction of 
multiple systems occurs very early on, within a few crossing 
times of the clusters (<lMyr, Parker et al. 2009) and, there¬ 
fore, at the age of the associations studied here this process has 
long ended. 

As mentioned above, we can also investigate the relation¬ 
ship between the projected physical separation and the age of 
the systems. Although the frequency of systems has not signifi¬ 
cantly evolved, it is possible that we can observe significant dif¬ 
ferences in the system parameters. To check this, we performed 
a KS test for three sets of ages: a) <30, b) 30-50, c) >50Myr. 
The age groups were comprised of 7, 23, 6 companions, respec¬ 
tively. The KS test produced one significant result; the distribu¬ 
tions of b) and c) are very unlikely to be resultant from two dif- 
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ferent parental distributions (KS statistic; 0.18, p-value; 0.99). 
This statistical similarity implies there is no significant dynami¬ 
cal evolution of the physical separations on this time-scale (30- 
100 Myr). Due to the insignificant results from the other KS tests 
we cannot discern whether or not the other age-grouped samples 
agree or not. 

6.4. Comparing physical separation distributions 

The current understanding of binary formation invokes two for¬ 
mation channels; (i) fragmentation of the proto-stellar core, sys¬ 
tems >100au (Bodenheimer & Burkert 2001) (ii) via gravita¬ 
tional instability and fragmentation of the proto-stellar accretion 
disc of the primary star, systems <100 au (Toomre 1964). If one 
of these two formation channels is dominant, one would expect 
to see a relative over- or under-abundance of systems between 
the two separation ranges, assuming that within our age range 
there is no significant evolution of the physical separations (see 
Section 6.3). 

Therefore, one interpretation of the over-abundance of mul¬ 
tiple systems with separations <100 au, a total of 21 systems, 
compared to wider systems, a total of 8, is that formation chan¬ 
nel (ii) is dominant for multiple systems in our sample. 

Another interpretation is that the original density of the 
SACY associations was much higher than the present day den¬ 
sity. This high-density environment would have dynamically 
processed many of the wider systems >100 au early-on and, 
therefore, would lead to an observed relative under-abundance in 
this separation range. However, this original density would have 
to be extremely high to account for such extreme processing. In 
addition to this, the MF value of SACY is large and similar to 
that of Taurus which is an indication that the binary population 
has not been heavily processed. 

6.5. Universality of multiple systems 

Using a compilation of multiplicity surveys. King et al. (2012a,b, 
K12a,b hereafter) claimed that binary (multiple system) forma¬ 
tion is not the same everywhere, and that star formation is not 
a single universal process. This conclusion was partly based on 
the significant differences found in the CDF between young pop¬ 
ulations and field stars, in the separation range of 19-100 au. 
They claim that, in this particular separation range, these sys¬ 
tems should not be dynamically processed and therefore their 
CDFs should be indistinguishable, independent of their age and 
initial environment. However, Marks et al. (2014) dispute many 
of the results and claims of K12b. They argue that if different 
initial densities are considered for the analysed SFRs, then the 
resultant multiple system properties are compatible with a sin¬ 
gle universal star formation process. In fact, the separation range 
19-lOOau is potentially not pristine, and in the case of Taurus, 
the comparable pristine range is only ^20-30 au. 

For the sake of comparison, we have also produced a sepa¬ 
ration distribution in the range 19-100 au, comparing 16 SACY 
systems, 23 Taurus systems, and 48 field systems. The results 
are shown in Figure 9. We note that we excluded data from D15, 
as in this separation range only seven systems were available for 
comparison. 

We find that the SACY and field distributions are statisti¬ 
cally indistinguishable within the CIs. We also performed a KS 
test which produced a p-value of 0.995, which is highly signif¬ 
icant, this means it is very unlikely that the two distributions 
are realisations of different parent distributions. If the separation 


range of 19-100 au is considered a pristine range, this result is a 
strong indication that the stars born in the field and those born in 
associations form in a similar way. 

In the case of Taurus, the distribution shows a different 
shape, and has a p-value of 0.30 and 0.24 when compared to 
SACY and the field, respectively. However, the p-values are not 
significant and therefore we cannot rule out that it is a realisation 
of the same parent distribution. 

SACY and Taurus share many similar multiplicity proper¬ 
ties (similar MFio-ioooau and TFio_ioooau values, a flat mass-ratio 
distribution). 

However, their respective separation distributions differ in 
both in the 10-1000 au and the 19-100 au range and additionally 
their MFio-ioooau,Mi>o .7 are marginally discrepant. If one pop¬ 
ulation had been more heavily dynamically processed than the 
other, we would expect to see differences in both their MF val¬ 
ues and their separation distributions. We do see tentative differ¬ 
ences in both quantities, implying that Taurus has experienced 
less dynamical processing than the SACY associations. 

At this time, the significance of the differences between the 
populations is still low (ail cr), and therefore further data and 
analysis is still needed to clarify whether these differences are 
physical or merely an effect of low-number statistics. One way 
to approach this problem is to focus on imaging a large num¬ 
ber of targets in a pristine separation range. However, calculat¬ 
ing this pristine range can be problematic due to unknown ini¬ 
tial densities. Once again, the recently launched ESA mission, 
GAIA, will put important constraints on the galactic velocities of 
members with accurate astrometry. This will provide much bet¬ 
ter constraints on the initial densities of such associations. 

7. Conclusions 

In this work we have presented multiplicity statistics of the 
SACY associations from AO-observations. We have derived 
multiplicity properties of our sample from observations of 113 
targets in the mass range 0.2-1.2 Mq. We compared these derived 
properties to other populations searching for statistical similari¬ 
ties and / or differences. 

This study and previous work (Elliott et al. 2014) form part 
of a larger project to characterise the multiplicity properties of 
the SACY associations across a large and continuous parame¬ 
ter space, using spectroscopy, sparse-aperture masking (SAM), 
speckle imaging, AO-imaging and classical direct imaging. 

Our observed sample consisted of 201 targets; 113 confirmed 
SACY members and 88 targets belonging to other populations 
(both from the field and PMS regions). As discussed in Section 3, 
these 88 targets were not included in our analysis. The parame¬ 
ters from individual observations can be found in Table B.l, and 
a summary of the identified multiple systems in Table C.l. 

Erom the multiplicity study of the 113 confirmed SACY 
members our main conclusions can be summarised as follows; 

1. We have identified 31 multiple systems (28 binaries, 3 
triples) from observations of 113 targets. 

2. Out of these 31 systems, 7 were identified from co-moving 
proper motion analysis, and 24 from contamination likeli¬ 
hood based on the component’s parameters (with probabili¬ 
ties >0.95). 

3. The multiplicity frequency and triple frequency for SACY, 
in the separation range 3-1000 au, are ME3_ioooau=28.4^3 g% 
andTE3_ioooau=2.8j;2-5 %, respectively. 

4. In the separation range 10-1000 au and for primary masses 
>0.7 Mo the multiplicity frequency (MEio_ioooau,Mi>o.7) of 
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Table 5. A summary of the properties of multiple systems identified in this work. 


ID 

Ang. sep. 
(") 

PA 

(deg) 

Dist. 

(pc) 

Phys. sep 
(au) 

AK 

(mag.) 

Comp. 

Prob.*"' 

Pop. 

Ml 

(Mo) 

Mi 

(Mo) 

q 

O-q 

GJ4231 

0.16 

298.4 

30.6 

4.8 

0.58 

A, B 

.99 

ABD 

0.5 

0.39 

0.78 

0.2 

GJ4231 

0.54 

0.4 

30.6 

16.6 

0.70 

A,C 

.99 

ABD 

0.5 

0.36 

0.73 

0.2 

TYC 91-82-1 

0.74 

110.4 

87.6 

64.8 

1.67 

A, B 

.99 

ABD 

0.82 

0.49 

0.6 

0.13 

BD-03 4778 

2.53 

208.5 

71.5 

180.9 

2.19 

A, B 

- 

ABD 

0.91 

0.46 

0.51 

0.12 

Wolf 1225 

1.44 

232.7 

15.1 

21.7 

0.07 

A, B 

- 

ABD 

0.48 

0.47 

0.99 

0.2 

HD 217379 

2.24 

239.0 

32.7 

73.2 

0.63 

A, B 

.99 

ABD 

0.81 

0.67 

0.82 

0.18 

HD 152555 

3.79 

58.0 

46.8 

177.4 

3.89 

A, B 

.77 

ABD 

1.13 

0.28 

0.25 

0.07 

GJ 3322 

1.36 

149.5 

37.8 

51.4 

0.78 

A, B 

.99 

BPC 

0.85 

0.57 

0.67 

0.18 

GSC 08350-01924 

0.73 

19.9 

66.3 

48.4 

0.17 

A, B 

- 

BPC 

0.63 

0.56 

0.9 

0.19 

2MAS S J05203182-t-0616115 

0.42 

235.0 

71.0 

29.8 

2.19 

A, B 

.96 

BPC 

0.66 

0.15 

0.23 

0.05 

CD-27 11535 

0.08 

278.8 

87.3 

7.0 

0.17 

A, B 

- 

BPC 

1.01 

0.97 

0.95 

0.2 

CD-26 13904 

0.27 

81.8 

78.9 

21.3 

0.93 

A, B 

- 

BPC 

1.03 

0.74 

0.72 

0.21 

2MASS J09131689-5529032 

0.14 

128.7 

131.7 

18.4 

1.48 

A, B 

.99 

CAR 

1.15 

0.72 

0.63 

0.09 

CD-57 1709 

5.84 

220.5 

100.5 

586.9 

2.28 

A, B 

.96 

CAR 

0.87 

0.34 

0.39 

0.15 

2MASS J08371096-5518105 

5.72 

23.5 

180.8 

1034.2 

3.79 

A, B 

.38 

CAR 

1.11 

0.18 

0.16 

0.04 

HD 22213 

1.67 

281.1 

52.4 

87.5 

1.87 

A, B 

.99 

CAR 

0.94 

0.54 

0.58 

0.14 

BD-07 2388 

0.11 

328.6 

42.6 

4.7 

0.62 

A, B 

.99 

CAR 

0.77 

0.64 

0.83 

0.21 

CPD-62 1293 

0.13 

344.9 

69.4 

9.0 

1.62 

A, B 

.98 

CAR 

0.7 

0.28 

0.41 

0.14 

V* V1221 Tau 

0.92 

134.3 

81.5 

74.6 

0.57 

A, B 

.99 

COL 

0.92 

0.81 

0.88 

0.15 

V* V1221 Tau 

2.17 

320.7 

81.5 

176.6 

2.5 

A,C 

.96 

COL 

0.92 

0.34 

0.37 

0.11 

2MASS J08240598-6334024 

0.17 

288.2 

118.6 

20.7 

3.141 

A, B 

.98 

COL 

1.2 

0.34 

0.29 

0.05 

2MASS J08240598-6334024 

6.38 

140.5 

118.6 

756.6 

3.358 

A,C 

.29 

COL 

1.2 

0.29 

0.25 

0.04 

2MASS J04272050-4420393 

0.43 

179.9 

86.2 

37.1 

2.15 

A, B 

.96 

COL 

0.81 

0.29 

0.36 

0.15 

2MASS J02303239-4342232 

0.13 

294.7 

51.2 

6.7 

1.09 

A, B 

.99 

COL 

0.81 

0.57 

0.71 

0.22 

BD-16 351 

0.34 

311.4 

80.9 

27.5 

0.66 

A, B 

.99 

COL 

0.83 

0.71 

0.85 

0.2 

CD-43 1395 

0.26 

324.6 

141.0 

36.7 

0.02 

A, B 

.99 

COL 

0.93 

0.92 

0.99 

0.16 

HD 272836 

1.52 

358.2 

78.5 

119.4 

2.492 

A, B 

.96 

COL 

0.84 

0.26 

0.31 

0.13 

HD 105923 

1.96 

145.1 

116.9 

229.1 

2.81 

A, B 

- 

ECH 

1.2 

0.39 

0.33 

0.02 

2MASS J12202177-7407393 

0.23 

3.7 

114.8 

26.4 

3.91 

A, B 

.48 

ECH 

1.07 

0.08 

0.07 

0.02 

TYC 9245-535-1 

0.11 

9.3 

119.3 

13.1 

1.06 

A, B 

.99 

ECH 

1.14 

0.63 

0.55 

0.1 

CD-30 3394A 

12.06 

116.7 

162.4 

1958.5 

0.5 

A, Ba 

.82 

OCT 

1.2 

1.14 

0.95 

0.08 

CD-30 3394A 

0.11 

43.2 

162.4 

17.9 

0.18 

Ba, Bb 

.99 

OCT 

1.14 

1.10 

0.96 

0.14 

HD 274576 

3.14 

15.8 

117.1 

367.7 

2.51 

A, B 

.95 

OCT 

0.89 

0.31 

0.35 

0.12 

2MASS J04302731-4248466 

0.49 

193.4 

122.8 

60.2 

0.13 

A, B 

.99 

OCT 

0.82 

0.8 

0.97 

0.2 

BD-20 1111 

4.01 

301.0 

129.7 

520.1 

3.56 

A, Ba 

.43 

OCT 

1.11 

0.21 

0.19 

0.05 

BD-20 1111® 

0.07 

340.0 

129.7 

9.1 

0.46 

Ba, Bb 

.98 

OCT 

0.21 

0.18 

0.86 

0.14 

CD-47 1999 

0.22 

321.5 

166.7 

36.7 

3.56 

A, B 

.81 

OCT 

1.2 

0.3 

0.25 

0.04 

2MASS J06033540-4911256 

0.75 

90.0 

173.8 

130.4 

0.06 

A, B 

.99 

OCT 

0.89 

0.88 

0.99 

0.18 

BD-18 4452 

3.18 

190.7 

16.7 

53.1 

0.81 

A, B 

- 

OCT 

0.51 

0.29 

0.58 

0.13 

BD-18 4452 

0.31 

349.4 

16.7 

5.2 

0.01 

Ba, Bb 

- 

OCT 

0.29 

0.28 

0.99 

0.2 

TYC 8083-455-1 

1.06 

347.8 

54.8 

58.1 

1.15 

A, B 

.98 

THA 

0.7 

0.35 

0.5 

0.14 

2MASS J05182904-3001321 

0.69 

76.5 

66.2 

45.7 

2.48 

A, B 

.97 

THA 

0.76 

0.17 

0.22 

0.07 

CD-44 1173 

0.42 

83.4 

44.0 

18.5 

2.44 

A, B 

.96 

THA 

0.75 

0.17 

0.22 

0.07 

HD 22705 

2.71 

89.2 

43.3 

117.3 

3.39 

A, B 

.52 

THA 

1.04 

0.22 

0.21 

0.06 


Notes. Probabilities are shown unless system was confirmed by co-moving proper motion analysis.**' Ba, Bb has a probability >0.95, however, 
currently there is no information regarding its membership and is therefore not included in further analysis. 


SACY (23.2+4-^%) is similar to the field (28.9+3 °%) and dis¬ 
crepant at a 1 cr level to Taurus (37.9+° ^%). 

5. The multiplicity frequency is not a function of the primary 
mass, within the primary mass range of 0.2-1.2Mo. 

6 . There is no evidence of dynamical destruction of multiple 
systems in the age range of lO-lOOMyr. 

7. There is evidence that the multiple systems’ physical sep¬ 
arations do not evolve significantly on the time-scale 30- 
100 Myr. 

8 . The mass-ratio distribution (power-law index 'y=-0.04 + 
0.14) is compatible with a uniform distribution (7=0). 

9. The SACY and the field separation distributions are indis¬ 
tinguishable in the separation range 19-100 au, suggesting a 


similar formation mechanism (if we assume this separation 
range is free of dynamical processing). 

A detailed analysis of the statistical properties of the SACY 
associations in the context of star formation will be presented 
in a future paper. It will include Elliott et al. (2014), this paper, 
and (Torres et al. 2008) in addition to collating all the available 
published information and data from the literature (e.g. Evans 
et al. 2012). We can then compare a significant number of these 
young systems across a larger parameter space, and make our 
analysis much more comprehensive. 
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Appendix A: All sources observed in this work 


Table A.l. All targets analysed in this work. 


ID 

RA 

hh:mm:ss.s 

DEC 

dd:mm:ss 

V 

(mag.) 

Ks 

(mag.) 

SpT 

(mas) 

1^6 

(mas) 

Plx 

(pc) 

Population 

Simbad Class. 

Lit. Flag" 

SACY Flag'' 

HD 1405 

00 18 20.9 

+30 57 22 

8.86 

6.51 

K2V 

143.7 

-171.5 

- 

ABD 

RS* 



HD 24681 

03:55:20.4 

-01:43:45 

9.05 

7.25 

G5 

46.2 

-90.2 

~ 

ABD 




TYC 91-82-1 

04:37:51.4 

+05:03:08 

11.01 

8.65 

- 

18.50 

-59.30 

~ 

ABD 



AB 

V* PX Vir 

13:03:49.6 

-05:09:42 

7.69 

5.51 

G5V 

-191.13 

-218.73 

46.10 

ABD 

BY* 


SBl 

LP 745-70 

16:33:41.6 

-09:33:11 

11.24 

7.55 

K9V 

-70.05 

-177.52 

32.60 

ABD 

PM* 

a*f' 


HD 160934 

17:38:39.6 

+61:14:16 

10.45 

6.81 

K7Ve 

-23.30 

47.71 

30.19 

ABD 

FI* 

AB® ,A*D, 













Aa,Ab' 


TYC 486-4943-1 

19:33:03.7 

+03:45:39 

11.15 

8.66 

- 

18.0 

-65.5 

~ 

ABD 




BD-03 4778 

20:04:49.3 

-02:39:20 

10.24 

7.92 

- 

28.6 

-71.7 

~ 

ABD 

Ro* 


AB 

HD 201919 

21:13:05.2 

-17:29:12 

10.64 

7.58 

K6Ve 

76.5 

-144.1 

~ 

ABD 

Ro* 



GJ4231 

21:52:10.4 

+05:37:35 

21.11 

7.38 

M2.4V 

119.17 

-150.29 

32.79 

ABD 

FI* 

AB,C 


GJ 885 A 

23:00:27.9 

-26:18:42 

10.48 

6.27 

(K7V)M+- 

113.2 

-169.2 

~ 

ABD 



SB3 

GJ 9809 

23:06:04.8 

+63:55:34 

10.87 

6.98 

M0.3V 

171.46 

-58.55 

40.81 

ABD 

FI* 

A*Y,Za*m' 


CD-45 14955B 

23:11:53.6 

-45:08:00 

13.8 

8.85 

M3Ve 

85.2 

-84.7 

19.7 

ABD 




HD 222575 

23:41:54.2 

-35:58:39 

9.39 

7.62 

G8V 

69.49 

-67.53 

15.70 

ABD 

Ro* 



HD 199058 

20:54:21.0 

+09:02:23 

8.62 

6.96 

G5 

37.2 

-56.9 

~ 

ABD'" 




TYC 1090-543-1 

20:54:28.0 

+09:06:06 

11.68 

8.82 

- 

34.4 

-58.6 

~ 

ABD'" 

V* 



HD 189285 

19:59:24.1 

-04:32:06 

9.43 

7.84 

G5 

14.2 

-55.8 

~ 

ABD'" 




V* PW And 

00:18:20.8 

+30:57:22 

8.81 

6.39 

K2V 

143.7 

-171.5 

~ 

ABD'* 

RS* 



V* LO Peg 

21:31:01.7 

+23:20:07 

9.19 

6.38 

K5-7V 

133.38 

-145.24 

40.32 

ABD'* 

BY* 

AB,AC' 


BD-i-41 4749 

23:19:39.5 

+42:15:09 

8.93 

7.22 

GO 

77.52 

-66.90 

19.94 

ABD'* 




HD 152555 

16:54:08.1 

-04:20:24 

7.82 

6.36 

F8/G0V 

-37.25 

-114.05 

21.40 

ABD 


AB 


Wolf 1225 

22:23:29.0 

+32:27:32 

11.45 

6.05 

MIV 

255.3 

-207.8 

- 

ABD, PMG'^ 


AB^ 

AB 

CD-52 9381 

20:07:23.7 

-51:47:27 

10.5 

7.39 

K6Ve 

88.9 

-146.2 

~ 

ARG 

Ro* 



GJ 3322 

05:01:58.7 

+09:58:59 

11.98 

6.37 

M3V 

12.09 

-74.41 

30.12 

BPC 

Ro* 

AB^ 

AB 

2MASS J05200029-I-0613036 

05:20:00.2 

+06:13:03 

11.58 

8.57 

K3 

5.5 

-37.8 

~ 

BPC 




2MASS J05203182+0616115 

05:20:31.8 

+06:16:11 

11.68 

8.57 

K3 

12.4 

-38.8 

~ 

BPC 



AB 

CD-27 11535 

17:15:03.6 

-27:49:39 

11.12 

7.38 

K5Ve 

-9.30 

-48.30 

~ 

BPC 

Ro* 


AB, SBl 

GSC 08350-01924 

17:29:20.6 

-50:14:53 

13.47 

7.99 

M3Ve 

-5.0 

-54.7 

~ 

BPC 



AB 

CD-54 7336 

17:29:55.0 

-54:15:48 

9.55 

7.36 

KIV 

-7.0 

-63.1 

~ 

BPC 

Ro* 



HD 161460 

17:48:33.7 

-53:06:43 

9.61 

6.78 

KOIV 

-5.60 

-53.30 

~ 

BPC 

Ro* 



V* V4046 Sgr 

18:14:10.4 

-32:47:34 

10.94 

7.25 

K5V 

2.1 

-54.5 

~ 

BPC 

BY* 


SB2 

GSC 07396-00759 

18:14:22.0 

-32:46:10 

12.78 

8.54 

MlVe 

7.3 

-39.9 

~ 

BPC 

Em* 



Smethells 20 

18:46:52.5 

-62:10:36 

12.08 

7.85 

MlVe 

13.7 

-82 

~ 

BPC 

pr* 



CD-31 16041 

18:50:44.4 

-31:47:47 

11.2 

7.46 

K8Ve 

10.6 

-77.8 

- 

BPC 

Ro* 

A*j' 


TYC 6872-1011-1 

18:58:04.1 

-29:53:04 

11.78 

8.02 

MOVe 

6.7 

-42.9 

~ 

BPC 

Ro* 



CD-26 13904 

19:11:44.6 

-26:04:08 

10.39 

7.37 

K4V(e) 

18.9 

-44.2 

- 

BPC 

pr* 

AB^ 

AB 

CPD-72 2713 

22:42:48.9 

-71:42:21 

10.6 

6.89 

K7Ve 

94.1 

-54.4 

~ 

BPC 

Ro* 



V* WW PsA 

22:44:57.9 

-33:15:01 

12.07 

6.93 

M4IVe 

184.76 

-119.76 

42.84 

BPC 

BY* 

A^ 


V* TX PsA 

22:45:00.0 

-33:15:25 

13.36 

7.79 

M5IVe 

183 

-118 

- 

BPC 

FI* 

B^ 


BD-13 6424 

23:32:30.8 

-12:15:51 

10.5 

6.57 

MOVe 

139.20 

-83.40 

~ 

BPC 

Ro* 



HD 8813 

01:23:25.8 

-76:36:42 

8.46 

6.79 

G6V 

100.16 

-20.93 

21.52 

CAR 

pr* 



HD 22213 

03:34:16.3 

-12:04:07 

8.86 

6.79 

G8V 

80.7 

-35.4 

~ 

CAR 

Ro* 


AB 

CD-44 1533 

04:22:45.6 

-44:32:51 

10.47 

8.58 

KOV 

32.0 

2.8 

~ 

CAR 

Ro* 



CD-57 1709 

07:21:23.7 

-57:20:37 

10.73 

8.7 

KOV 

-1.6 

22.1 

~ 

CAR 

Ro* 


AB 

TYC 8557-1251-1 

07:55:31.6 

-54:36:50 

11.44 

9.19 

G9V 

-4.1 

10.3 

~ 

CAR 

Ro* 



2MASS J08110934-5555563 

08:11:09.3 

-55:55:56 

11.52 

9.4 

G8V(e) 

-3.5 

15.1 

~ 

CAR 

pr* 




P. Elliott et al.: Search for associations containing young stars (SACY), Online Material p 1 



ID 


RA DEC V mag. Kj mag. SpT jis Fix Population Simbad Class. Lit. Flag" SACY Flag'* 


BD-07 2388 

08:13:50.9 

-07:38:24 

9.32 

6.92 

KO 

-25.3 

-45.2 

~ 

CAR 

Ro* 

AB 

2MASS J08371096-5518105 

08:37:10.9 

-55:18:10 

11.04 

9.38 

G9V 

-7.4 

18.0 

~ 

CAR 

pr* 

AB? 

CD-61 2010 

08:42:00.4 

-62:18:26 

10.95 

8.83 

KOV 

-11.3 

15.9 

~ 

CAR 

Ro* 


CD-53 2515 

08:51:56.4 

-53:55:56 

11.06 

8.75 

G9V 

-11.6 

12.7 

~ 

CAR 

Ro* 


2MASS J08521921-6004443 

08:52:19.2 

-60:04:44 

11.48 

9.37 

KOV 

-2.8 

12.0 

~ 

CAR 

pr* 


2MASS J08563149-5700406 

08:56:31.4 

-57:00:40 

11.83 

9.76 

KOV(e) 

-14.7 

6.0 

~ 

CAR 

pr* 


TYC 8582-3040-1 

08:57:45.6 

-54:08:36 

11.71 

9.35 

K2IV(e) 

-13.1 

8.6 

~ 

CAR 

Ro* 


CD-49 4008 

08:57:52.1 

-49:41:50 

10.5 

8.64 

G9V 

-23.00 

15.40 

~ 

CAR 

Ro* 


CD-54 2499 

08:59:28.7 

-54:46:49 

10.08 

8.4 

G5IV 

-18.3 

15.2 

~ 

CAR 

Ro* 


CD-55 2543 

09:09:29.3 

-55:38:27 

10.21 

8.4 

G8V 

-13.2 

14.8 

~ 

CAR 

Ro* 


TYC 8174-1586-1 

09:11:15.8 

-50:14:14 

11.81 

9.5 

K5Ve 

-19.45 

7.42 

~ 

CAR 

Ro* 


2MASS J09131689-5529032 

09:13:16.8 

-55:29:03 

11.36 

8.36 

G5V(e) 

-13.8 

14.8 

~ 

CAR 

pr* 

AB 

HD 309751 

09:31:44.7 

-65:14:52 

11.33 

8.36 

G3V 

-12.40 

19.00 

~ 

CAR 

pr* 

SB2 

CPD-62 1293 

09:43:08.8 

-63:13:04 

10.44 

8.6 

G6V 

-34.9 

25.0 

~ 

CAR 

Ro* 

AB 

TYC 9217-417-1 

09:59:57.6 

-72:21:47 

11.7 

8.69 

K4Ve 

-24.2 

29.0 

~ 

CAR 

Ro* 


CD-69 783 

10:41:23.0 

-69:40:43 

10.27 

8.38 

G8V(e) 

-30.0 

24.4 

~ 

CAR 

Ro* 


TYC 8962-1747-1 

11:08:07.9 

-63:41:47 

12.07 

8.29 

MOVe 

-33.2 

7.2 

~ 

CAR 

Ro* 


HD 107722 

12:23:29.0 

-77:40:51 

8.24 

7.14 

F6V 

-65.1 

10.8 

~ 

CAR 

-ppH; 


TYC 9486-927-1 

21:25:27.4 

-81:38:27 

11.85 

7.34 

MlVe 

63.4 

-110.0 

~ 

CAR 

Ro* 


BD-03 5579 

23:09:37.1 

-02:25:55 

10.99 

7.82 

K4Ve 

57.5 

-47.4 

~ 

CAR 

Ro* 


CPD-53 1875 

08:45:52.7 

-53:27:28 

10.42 

8.77 

G2V 

-13.8 

7.7 

~ 

car'“ 

Ro* 

AB 

CPD-55 1885 

09:00:03.3 

-55:38:24 

10.83 

9.01 

G5V 

-12.7 

16.4 

~ 

car'“ 

Ro* 


Ass ChaT 2-21 

11:06:15.4 

-77:21:56 

11.42 

6.42 

G5Ve 

-16.5 

-0.5cm 

~ 

CHA I 

pr* 

D+w® 

2MASS J12571172-7640111 

12:57:11.7 

-76:40:11 

15.0 

8.4 

MOe 

10 

-32 

- 

CHA II'‘ 

Y*0 

AB 

BD-16 351 

02:01:35.6 

-16:10:00 

10.34 

7.96 

- 

57.1 

-29.7 

~ 

COL 

Ro* 

AB 

2MASS J02303239-4342232 

02:30:32.4 

-43:42:23 

10.37 

7.23 

K5V(e) 

80.5 

-14.9 

~ 

COL 

pr* 

AB 

BD-11 648 

03:21:49.6 

-10:52:17 

11.63 

9.26 

G9 

~ 

~ 

~ 

COL 

Ro* 


V* V1221 Tau 

03:28:14.9 

+04:09:47 

9.7 

7.44 

G6V+K3V... 

34.4 

-45.0 

~ 

COL 

BY* 

AB,C 

2MASS J03573723-0416159 

03:57:37.2 

-04:16:15 

10.72 

8.75 


~ 

~ 

~ 

COL 

Ro* 


HD 26980 

04:14:22.5 

-38:19:01 

9.08 

7.62 

G3V 

41.78 

1.40 

12.42 

COL 

Ro* 


HD 27679 

04:21:10.3 

-24:32:20 

9.43 

7.81 

G2V 

40.0 

-9.1 

~ 

COL 

pr* 


CD-43 1395 

04:21:48.6 

-43:17:32 

10.18 

8.42 

G7V 

20.4 

5.6 

~ 

COL 

Ro* 

AB 

2MASS J04272050-4420393 

04:27:20.4 

-44:20:39 

10.91 

8.56 

KlV(e) 

31.3 

2.2 

~ 

COL 

pr* 

AB 

CD-36 1785 

04:34:50.7 

-35:47:21 

10.84 

8.59 

KlVe 

34.8 

-1.6 

~ 

COL 

Ro* 


BD-i-08 742 

04:42:32.0 

+09:06:00 

11.19 

9.12 

G5V: 

27.60 

-25.50 

~ 

COL 



HD 31242 

04:51:53.5 

-46:47:13 

9.85 

8.16 

G5V 

32.7 

14.2 

~ 

COL 

Ro* 

SBl 

HD 272836 

04:53:05.1 

-48:44:38 

10.72 

8.24 

K2V(e) 

30.8 

17.0 

~ 

COL 

Ro* 


TYC 9178-284-1 

06:55:25.1 

-68:06:21 

11.91 

8.94 

K4Ve 

2.7 

20.2 

~ 

COL 

Ro* 


2MASS J08240598-6334024 

08:24:05.9 

-63:34:02 

9.87 

8.13 

G5V 

-15.6 

22.4 

~ 

COL 

pr* 

AB,C 

V* V479 Car 

09:23:34.9 

-61:11:35 

10.16 

7.96 

KlV(e) 

-28.30 

17.98 

11.36 

COL 

BY* 


2MASS J09322609-5237396 

09:32:26.0 

-52:37:39 

10.86 

8.84 

G8V(e) 

-16.5 

10.9 

~ 

COL 

pr* 


HD 174656 

18:53:05.9 

-36:10:22 

9.87 

7.28 

KOIV(e) 

0.20 

-21.40 

- 

CRA 

Ro* 

AB^ 

V* DZ Cha 

11:49:31.8 

-78:51:01 

12.0 

8.49 

MOVe 

-38.0 

-8.3 

~ 

ECH 

Or* 


2MASS J11594226-7601260 

11:59:42.2 

-76:01:26 

11.31 

8.3 

K4Ve 

-40.5 

-5.83 

9.89 

ECH 

-ppH; 


HD 104467 

12:01:39.1 

-78:59:16 

8.56 

6.85 

G3V(e) 

-42.8 

-4.0 

~ 

ECH 


SBl 

HD 105923 

12:11:38.1 

-71:10:36 

9.16 

7.18 

G8V 

-38.9 

-8.3 

~ 

ECH 

Ro* 

AB 

2MASS J12194369-7403572 

12:19:43.6 

-74:03:57 

12.5 

8.86 

Ml 

-42.1 

-11.7 

~ 

ECH 

-ppH: 


2MASS J12202177-7407393 

12:20:21.7 

-74:07:39 

12.0 

8.37 

MIV 

-41.0 

-6.5 

~ 

ECH 


AB? 

2MASS J12392124-7502391 

12:39:21.2 

-75:02:39 

10.3 

7.78 

K3Ve 

-43.6 

-11.1 

~ 

ECH 

-ppH: 


TYC 9245-535-1 

12:56:08.3 

-69:26:53 

12.06 

7.99 

K7Ve 

-33.6 

-7.6 

~ 

ECH 

Ro* 

AB 

CD-69 1055 

12:58:25.5 

-70:28:49 

9.67 

7.55 

KOVe 

-43.1 

-18.7 

~ 

ECH 

Ro* 


V* MP Mus 

13:22:07.5 

-69:38:12 

9.75 

7.29 

KlVe 

-40.80 

-23.30 

~ 

ECH 

-ppH: 


2MASS J12375231-5200055 

12:37:52.2 

-52:00:05 

10.68 

6.02 

M3V(e) 

-1032.41 

30.39 

103.18 

HYA‘2 

pr* 


2MASS J22420163+0946091 

22:42:01.6 

+09:46:09 

11.99 

8.26 

K6Ve 

17.27 

30.22 

26.73 

HYA‘2 

pr* 
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ID 

RA 

DEC 

V mag. 

Ks mag. 

SpT 

f^a 

^6 

Plx 

Population 

Simbad Class. 

Lit. Flag" 

SACY Flag'' 

2MASS J23350028+0136193 

23:35:00.2 

+01:36:19 

9.59 

6.04 

K7V 

341.31 

25.50 

52.56 

IC2391‘2 


A*F,BC' 


TYC 8639-1114-1 

11:55:42.9 

-56:37:31 

12.11 

8.04 

MOVe 

-45.4 

-6.0 

~ 

LCC 

Ro* 


AB 

TYC 8978-3494-1 

12:12:48.8 

-62:30:31 

12.13 

7.96 

K7Ve 

-36.4 

-16.4 

~ 

LCC 

Ro* 


AB 

TYC 8637-1558-1 

12:16:01.1 

-56:14:07 

11.22 

7.96 

K5Ve 

-42.7 

-10.3 

~ 

LCC 

Ro* 



CD-62 657 

12:28:25.3 

-63:20:58 

9.08 

7.33 

G7V 

-34.9 

-11.3 

~ 

LCC 

Ro* 



CPD-63 2367 

12:36:38.9 

-63:44:43 

9.69 

7.37 

K2V 

-37.8 

-9.7 

- 

LCC 

Ro* 

A*pl 


2MASS J12474824-5431308 

12:47:48.2 

-54:31:30 

11.82 

8.05 

MOVe 

-41.5 

-11.0 

~ 

LCC 

pr* 



HD 112245 

12:56:09.4 

-61:27:25 

9.63 

7.26 

KOVe 

-52.5 

-14.7 

~ 

LCC 

Ro* 


AB 

TYC 9012-1005-1 

13:44:42.7 

-63:47:49 

11.04 

7.74 

K4Ve 

-44.9 

-20.8 

~ 

LCC 

Ro* 

A,Z*r,A*Y' 


CPD-64 1859 

12:19:21.6 

-64:54:10 

9.87 

7.4 

K3V 

-40.20 

-10.20 

- 

LCC‘5 

Ro* 



TYC 7846-833-1 

15:56:44.0 

-42:42:29 

11.88 

8.27 

K6IVe 

-13.7 

-16.6 

~ 

LUP 

Ro* 


AB, C 

CD-33 11275 

16:35:22.4 

-33:28:53 

11.42 

8.08 

K4Ve 

-16.8 

-33.8 

~ 

LUP 

Ro* 


AB, C 

CD-58 860 

04:11:55.6 

-58:01:47 

10.01 

8.36 

G6V 

-7.9 

38.2 

~ 

OCT 

Ro* 



2MASS J04302731-4248466 

04:30:27.3 

-42:48:46 

10.75 

8.73 

G9V(e) 

2.6 

20.8 

~ 

OCT 

pr* 


AB 

HD 271037 

05:06:50.5 

-72:21:11 

10.91 

8.67 

KOIV 

-7.0 

24.7 

~ 

OCT 

Ro* 



HD 274576 

05:28:51.3 

-46:28:19 

10.5 

8.81 

G6V 

-10.3 

22.6 

~ 

OCT 

Ro* 


AB 

BD-20 1111 

05:32:29.3 

-20:43:33 

10.4 

8.7 

- 

-8.3 

12.0 

~ 

OCT 



Ba, Bb ? 

CD-47 1999 

05:43:32.1 

-47:41:10 

10.19 

8.64 

GOV 

-10.0 

15.0 

~ 

OCT 

Ro* 


AB? 

2MASS J05581182-3500496 

05:58:11.8 

-35:00:49 

11.24 

9.38 

G9V 

-12.7 

13.0 

~ 

OCT 

pr* 



2MASS J06033540-4911256 

06:03:35.4 

-49:11:25 

11.18 

9.08 

KOVe 

-9.3 

8.3 

~ 

OCT 

pr* 


AB 

CD-66 395 

06:25:12.3 

-66:29:10 

10.87 

9.0 

KOIV 

-18.40 

24.40 

~ 

OCT 

Ro* 



2MASS J06400573-3033089 

06:40:05.7 

-30:33:08 

10.24 

8.7 

F9V 

-12.60 

7.70 

~ 

OCT 

pr* 

AB^ 


BD-18 4452 

17:13:11.6 

-18:34:25 

10.86 

6.48 

MOVe 

39.0 

-88.7 

- 

OCT 

pr* 

ab' 

A, Ba, Bb 

CD-30 3394 

06:40:04.9 

-30:33:03 

9.83 

8.59 

F6V 

-18.3 

7.0 

~ 

OCT 



A, Ba, Bb 

V* BO Mic 

20:47:45.0 

-36:35:40 

9.39 

6.79 

K3V(e) 

11.42 

-75.87 

19.16 

PMG‘2 

BY* 



V* VI002 Sco 

16:12:40.5 

-18:59:28 

9.81 

7.49 

K6e 

-11.70 

-22.20 

~ 

RHO 

Or* 

AB,Aa' 

AB 

2MASS J16141107-2305362 

16:14:11.0 

-23:05:36 

10.5 

7.46 

K2IV 

-11.9 

-22.9 

~ 

RHO 


A*G,Aa,Ab' 


HD 147808 

16:24:51.3 

-22:39:32 

9.5 

7.08 

G9IVe 

-13.4 

-20.0 

- 

RHO 

-p-pH; 

AB,Aa,Ab' 


V* V896 Sco 

16:25:38.4 

-26:13:54 

11.75 

7.52 

M 

~ 

- 

~ 

RHO 

Or* 


AB 

2MASS J16312019-2430009 

16:31:20.1 

-24:30:00 

10.69 

7.09 

K3e 

-5.1 

-28.8 

~ 

RHO 

SB* 

AB 

AB, C 

V* V877 Sco 

16:33:41.9 

-25:23:34 

10.66 

7.95 

K4IVe 

-2.4 

-25.4 

~ 

RHO 

bL* 



2MASS J16471358-1514275 

16:47:13.5 

-15:14:27 

12.02 

8.04 

MOVe 

-5 

-11 

~ 

RHO'^ 

-p-pHi 


AB 

2MASS J16482187-1410427 

16:48:21.8 

-14:10:42 

14.79 

7.86 

- 

- 

- 

- 

RHO'^ 

Y*0 


AB 

CD-78 24 

00:42:20.3 

-77:47:39 

10.21 

7.53 

K3Ve 

79.5 

-32.1 

~ 

THA 

Ro* 



CD-53 544 

02:41:46.8 

-52:59:52 

10.22 

6.76 

K6Ve 

97.5 

-13.7 

~ 

THA 

Ro* 



CD-35 1167 

03:19:08.6 

-35:07:00 

11.0 

7.72 

K7Ve 

89.2 

-20.3 

~ 

THA 

Ro* 



CD-46 1064 

03:30:49.0 

-45:55:57 

9.79 

7.1 

K3V 

87.2 

-5.0 

~ 

THA 

Ro* 



CD-44 1173 

03:31:55.6 

-43:59:13 

10.9 

7.47 

K6Ve 

90.9 

-5.0 

~ 

THA 

Ro* 


AB 

HD 22705 

03:36:53.4 

-49:57:28 

7.65 

6.14 

G2V 

89.74 

0.29 

23.07 

THA 


AB,AC' 

AB 

BD-12 943 

04:36:47.1 

-12:09:20 

9.86 

7.76 

- 

45.9 

-20.4 

~ 

THA 

Ro* 



TYC 8083-455-1 

04:48:00.6 

-50:41:25 

11.53 

7.92 

K7Ve 

64.3 

16.6 

~ 

THA 

Ro* 


AB 

BD-19 1062 

04:59:32.0 

-19:17:41 

10.61 

8.07 

- 

42.6 

-10.0 

~ 

THA 

Ro* 



ASAS 1051536-0930.8 

05:15:36.4 

-09:30:51 

9.82 

8.08 

G5 

31.4 

-20.4 

~ 

THA 

-p-pH; 



2MASS J05182904-3001321 

05:18:29.0 

-30:01:32 

11.66 

8.3 

K4Ve 

37.7 

4.3 

~ 

THA 

pr* 


AB 

TYC 8098-414-1 

05:33:25.5 

-51:17:13 

11.74 

8.16 

K7Ve 

54.1 

16.6 

~ 

THA 

Ro* 



2MASS 105490656-2733556 

05:49:06.5 

-27:33:55 

10.22 

8.25 

G9V 

30.90 

-4.70 

~ 

THA 

pr* 



TYC 9344-293-1 

23:26:10.7 

-73:23:49 

11.83 

7.94 

MOVe 

71.5 

-66.5 

~ 

THA 

Ro* 



HD 30051 

04:43:17.2 

-23:37:42 

7.12 

6.02 

F2/3IV/V 

50.25 

-11.84 

15.73 

THA'® 




CD-46 10045 

15:25:11.6 

-46:59:12 

10.08 

7.78 

G9V 

-49.63 

-21.24 

15.46 

UCL 


AB^ 

AB 

V* KW Lup 

15:45:47.6 

-30:20:55 

9.28 

6.46 

K2V 

-69.11 

-98.40 

24.78 

ucl“> 

BY* 


AB 

2MASS 115594951-3628279 

15:59:49.5 

-36:28:27 

10.36 

8.03 

K3e 

-16.21 

-51.01 

15.40 

US 

pr* 



V* V857 Ara 

17:18:14.6 

-60:27:27 

9.48 

7.53 

G8V 

-54.62 

-91.04 

16.97 

US, THA"’ 

BY* 


AB 

TYC 6234-1287-1 

17:26:56.5 

-16:31:34 

11.59 

7.83 

K4Ve 

-10.1 

-39.0 

~ 

US?'^ 

Ro* 


AB, SB2? 
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ID 

RA 

DEC 

V mag. 

Ks mag. 

SpT 

f^a 

^6 

Plx 

Population 

Simbad Class. 

Lit. Flag" 

SACY Flag'' 

2MASS J09551508-6203324 

09:55:15.0 

-62:03:32 

12.01 

9.66 

K2V 

-15.3 

4.1 

~ 

N 

pr* 



BD+01 3657 

18:22:17.2 

+01:42:25 

10.2 

6.86 

K6Vke 

84.65 

-18.29 

37.28 

N 

Ro* 



HD 83096 

09:31:24.8 

-73:44:49 

7.3 

6.17 

FIV 

-35.40 

32.44 

11.98 

N 




HD 154361 

17:05:08.4 

-01:47:10 

9.5 

6.81 

KO 

75.7 

-44.2 

~ 

N 

cC* 



TYC 438-902-1 

18:03:17.8 

+04:48:26 

10.3 

8.01 

K2 

-11.8 

-27.0 

~ 

N 

Ro* 



TYC 458-290-1 

18:31:59.1 

+05:40:19 

10.82 

8.85 

- 

31.8 

-5.9 

- 

N 

Ro** 



TYC 460-624-1 

18:45:10.2 

+06:20:15 

10.76 

6.81 

- 

-38.0 

-83.7 

~ 

N 

Ro* 



TYC 1026-1952-1 

18:47:25.5 

+08:41:07 

12.1 

8.97 

- 

-11.1 

-18.2 

~ 

N 

Ro* 


AB 

BD+01 3828 

18:57:19.4 

+01:20:33 

9.56 

7.79 

GO 

27.4 

-0.9 

- 

N 


AB^ 


HD 176898 

19:02:17.0 

+02:44:21 

9.19 

7.16 

G 

24.6 

-30.2 

~ 

N 




TYC 490-110-1 

19:31:34.9 

+05:56:07 

11.13 

7.73 

K3III 

22.2 

2.8 

~ 

N 

Ro* 


AB 

TYC 482-1106-1 

19:36:04.3 

+03:25:48 

10.43 

8.33 

- 

3.0 

-24.7 

~ 

N 

Ro* 



TYC 1062-1904-1 

19:47:39.1 

+09:29:36 

11.81 

8.85 

- 

26.3 

13.2 

~ 

N 




TYC 1058-1925-1 

19:51:13.3 

+08:42:19 

10.48 

8.13 

- 

11.1 

-17.5 

~ 

N 



AB 

TYC 5188-193-1 

20:35:27.1 

-07:06:54 

11.59 

10.78 

- 

16.2 

-14.6 

~ 

N 




HD 199602 

20:58:32.5 

-08:50:22 

9.19 

6.4 

KO 

-31.9 

-56.1 

~ 

N 




TYC 538-573-1 

21:01:24.5 

+05:42:12 

11.55 

9.39 

- 

18.7 

-3.3 

~ 

N 

V* 


AB 

TYC 553-84-1 

21:56:22.0 

+05:17:48 

10.68 

6.9 

- 

-8.2 

-6.3 

~ 

N 




CCDM J00002+0146AB 

00:00:12.2 

+01:46:17 

10.78 

7.93 

G9IV 

44.10 

-22.90 

- 

N 


A2 

AB, SB3 

TYC 1-1187-1 

00:09:21.7 

+00:38:06 

11.88 

8.71 

K4Ve 

119.4 

-29.7 

~ 

N 

Ro* 


SB2 

HD 531 

00:09:51.6 

+08:27:11 

9.35 

7.57 

G6V+G7V 

56.21 

-11.41 

13.11 

N 


AB^ 

AB 

GJ 3029 

00:21:37.2 

-46:05:33 

12.22 

7.45 

M3Ve 

-302.66 

-360.97 

51.80 

N 

FI* 



BD+07 85 

00:38:59.2 

+08:28:41 

10.09 

7.2 

K4Ve 

100.9 

-48.9 

~ 

N 

Ro* 



2MASS J01452133-3957204 

01:45:21.3 

-39:57:20 

11.39 

7.59 

MIV 

279.02 

131.42 

33.84 

N 

pr* 

ab’ 


BD-21 332 

01:53:11.3 

-21:05:43 

11.73 

7.14 

M2Ve 

271.4 

72.5 

- 

N 

Ro* 

SB2* 


GJ3148 A 

02:16:41.1 

-30:59:18 

12.05 

7.13 

M3Ve 

683.07 

250.88 

69.81 

N 

FI* 

AB^ 


2MASS J14110874-6155469 

14:11:08.7 

-61:55:47 

10.8 

6.92 

MlVe 

147.9 

-51.0 

~ 

N 

pr* 



2MASS J15011648-4339311 

15:01:16.5 

-43:39:31 

10.8 

6.07 

M3Ve 

98.0 

26.9 

~ 

N 

pr* 


SB2 

2MASS J15104047-5248189 

15:10:40.4 

-52:48:19 

12.26 

6.85 

M2V(e) 

-145.3 

-167.5 

~ 

N 

pr* 


AB, SB3 

2MASS J15240033-7054096 

15:24:00.3 

-70:54:09 

11.93 

7.77 

MlVe 

-145.1 

-128.2 

~ 

N 

pr* 


SB2 

CD-24 12794 

16:43:56.9 

-25:08:36 

11.26 

7.42 

K6Ve 

8.7 

-50.7 

~ 

N 

Ro* 



TYC 6242-104-1 

17:21:56.0 

-20:10:51 

13.22 

9.15 

K5Ve 

-11.7 

-13.7 

~ 

N 

Ro* 



V* V2384 Oph 

17:42:30.3 

-28:44:55 

9.03 

7.32 

G5V 

-2.83 

-40.87 

12.82 

N 

bL* 



HIP 96515 

19:37:08.7 

-51:34:00 

12.29 

7.96 

MlVe 

91.07 

-21.94 

22.79 

N 

Al* 

AB^ 

AB, SB2 

GJ 781.1 A 

20:07:44.9 

-31:45:14 

13.08 

7.4 

M4Ve 

286.55 

-750.5 

67.08 

N 

FI* 

AB^ 

SB2 

HIP 106043 

21:28:44.4 

-47:15:42 

12.2 

8.29 

MlVe 

-16.42 

-58.10 

27.94 

N 

Ro* 



GJ 841 A 

21:57:41.2 

-51:00:22 

10.5 

5.88 

M2Ve 

-35.66 

-377.25 

62.61 

N 

FI* 

AB^ 


BD+02 4456 

21:59:59.9 

+03:02:24 

10.35 

7.89 

KOIV 

63.80 

-3.20 

~ 

N 

Ro* 


SB2 

2MASS J22023015-0406117 

22:02:30.1 

-04:06:11 

10.06 

7.97 

K2V 

29.9 

-13.4 

~ 

N 

pr* 



V* CS Gru 

22:15:35.2 

-39:00:50 

9.32 

7.09 

KOV 

94.36 

-39.79 

17.84 

N 

BY* 


SBl 

HD 212781 

22:26:53.2 

-00:50:39 

10.46 

8.34 

KOV(e) 

-4.30 

-48.00 

~ 

N 

Ro* 


AB, SBl 

HIP 113001 

22:53:05.9 

+07:28:19 

11.72 

7.94 

K6Ve 

59.38 

-79.98 

19.62 

N 

Ro* 


SB2 

TYC 572-382-1 

22:56:49.5 

+02:35:39 

10.29 

8.08 

KlIV 

-14.2 

-18.0 

~ 

N 

Ro* 


AB 

2MASS J23045301+0632003 

23:04:53.0 

+06:32:00 

10.99 

8.18 

KlVe 

-14.8 

-23.4 

~ 

N 

pr* 



2MASS J23085046+0000528 

23:08:50.4 

+00:00:52 

10.5 

8.62 

G8V 

16.8 

-7.1 

~ 

N 

pr* 



TYC 584-343-1 

23:21:56.3 

+07:21:32 

10.97 

8.94 

KOV 

21.50 

-2.50 

~ 

N 

Ro* 


AB 

V* BS Psc 

23:52:24.3 

-00:55:59 

10.76 

8.95 

G7V 

19.90 

-12.30 

~ 

N 

BY* 



CD-36 14261 

20:36:08.3 

-36:07:11 

11.73 

7.17 

M3Ve 

11.4 

40.3 

~ 

N 

Ro* 
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(2014); (15) Mamajek et al. (2006); (16) Song et al. (2012); (17) Sartori et al. (2003); (18) Lopez-Santiago et al. (2006); (19) Zuckerman et al. (201 1); (20) Zuckerman et al. (2004); (21) Schlieder 
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Notes. Formats as ’ A*Y’ signifies multiple components AB, AC, ..., AY Note in SIMBAD to potential component within 30" Note in SIMBAD for triple system identified by Malogolovets 
et al. (2007) SB systems denote spectroscopic binary systems identified in Torres et al. (2008); Elliott et al. (2014). 
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Appendix B: Individual detections in this work 

Table B.l. Details of individual detections in this work. 
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Table C.l. A summary of the properties of non-SACY multiple systems identified in this work. 


ID 

Ang. sep. 
(") 

PA 

(deg.) 

Dist. 

(pc) 

Phys. sep 
(au) 

AK 

(mag.) 

Comp. 

Prob. 

Pop. 

q 

O-q 

PMS populations 

CPD-53 1875 

0.21 

288.6 

_ 

_ 

2.18 

A, B 

0.91 

CAR* 

0.27 

0.08 

2MASS J12571172-7640111 

0.28 

150.1 

- 

- 

0.69 

A, B 

- 

CHAlP 

0.61 

0.05 

HD 112245 

2.56 

169.4 

66.0 

168.9 

2.48 

A, B 

- 

LCA, C 

0.22 

0.12 

TYC 8978-3494-1 

0.43 

41.2 

72.5 

31.2 

0.13 

A, B 

- 

LCA, C 

0.93 

0.2 

TYC 8639-1114-1 

0.84 

305.6 

79.0 

66.4 

0.42 

A, B 

- 

LCA, C 

0.80 

0.2 

CD-33 11275 

0.18 

252.9 

123.7 

22.3 

0.04 

A, B 

- 

LUP 

1.03 

0.2 

CD-33 11275 

1.78 

3.7 

123.7 

220.2 

0.54 

A, C 

- 

LUP 

0.67 

0.13 

TYC 7846-833-1 

0.2 

192.5 

139.7 

27.9 

0.76 

A, B 

- 

LUP 

0.58 

0.11 

TYC 7846-833-1 

6.0 

19.8 

139.7 

837.9 

2.47 

A,C 

- 

LUP 

0.17 

0.04 

V* V896 Sco 

0.87 

20.9 

125.0 

108.7 

1.68 

A, B 

- 

RHO 

0.30 

0.05 

V* V1002 Sco 

0.12 

156.0 

138.7 

16.6 

1.58 

A, B 

0.89 

RHO 

0.32 

0.03 

2MASS J16312019-2430009 

0.32 

170.4 

103.4 

33.1 

2.98 

A, B 

- 

RHO 

0.12 

0.05 

2MASS J16312019-2430009 

4.48 

12.1 

103.4 

463.2 

0.33 

A,C 

- 

RHO 

0.79 

0.14 

2MASS J16471358-1514275 

0.35 

60.5 

- 

- 

1.15 

A, B 

- 

RHO^ 

0.44 

0.05 

2MASS J16482187-1410427 

1.01 

252.2 

- 

- 

1.93 

A, B 

- 

RHO* 

0.26 

0.05 

CD-46 10045 

0.62 

334.5 

67.9 

42.1 

0.07 

A, B 

- 

UCL 

0.96 

0.2 

V* V857 Ara 

0.2 

194.2 

- 

- 

2.45 

A, B 

- 

US 

0.19 

0.06 

TYC 6234-1287-1 

1.51 

245.2 

- 

- 

0.21 

A, B 

- 

US?^ 

0.88 

0.05 

Other 

HIP 96515 

8.58 

224.0 

_ 

_ 

6.36 

A, B 

_ 

N, Lield 

_ 

_ 

HD 212781 

5.34 

77.3 

- 

- 

4.71 

A, B 

- 

N, Lield 

- 

- 

TYC 1026-1952-1 

6.99 

58.6 

- 

- 

3.76 

A, B 

- 

N, Lield 

- 

- 

TYC 572-382-1 

3.07 

53.5 

- 

- 

2.68 

A, B 

- 

N, Lield 

- 

- 

TYC 490-110-1 

2.42 

180.3 

- 

- 

2.45 

A, B 

- 

N, Lield 

- 

- 

TYC 538-573-1 

6.19 

25.3 

- 

- 

2.27 

A, B 

- 

N, Lield 

- 

- 

TYC 584-343-1 

5.08 

228.5 

- 

- 

1.53 

A, B 

- 

N, Lield 

- 

- 

2MASS 115104047-5248189 

0.47 

294.7 

- 

- 

1.14 

A, B 

- 

N, Lield 

- 

- 

CCDM 109314-7345AB 

0.95 

133.1 

- 

- 

0.62 

A, B 

0.99 

N, Lield 

- 

- 

TYC 1058-1925-1 

0.76 

236.1 

- 

- 

0.17 

A, B 

- 

N, Lield 

- 

- 

CCDM 100002+0146AB 

1.81 

262.9 

- 

- 

0.16 

A, B 

- 

N, Lield 

- 

- 

HD 531 

5.21 

275.9 

- 

- 

0.05 

A, B 

- 

N, Lield 

- 

- 


References. (1) Messina et al. (2010); (2) Alcala et al. (2008); (3) Sartori et al. (2003); (4) Aarnio et al. (2008) 
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